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PREFACE 

This report presents the results of an experiment in remote 
bathymetry conducted jointly as a NASA/Cousteau Society Ocean Bathymetry 
experiment in August and September of 1975. Fabian C. Polcyn, Senior 
Research Engineer of the Environmental Research Institute of Michigan, 
served as Principal Investigator. Many personnel from several organi- 
zations participated in the preparation and execution of the experiment. 
Among these are: 

^ Cousteau Society 

Captain Jacques Ccusteau, Experiment Co-Sponsor 
* NASA Headquart ers 

Mr. Russell Schweickart, Experiment Coordinator 
^ NASA Goddard 

Dr. Enrico P. Mercanti, Experiment Manager 
Mr. Charles Bohn, Experiment Manager and Technical 
Monitor 

Dr. Ross McCluney, Experiment Scientific Monitor 
Mr. Charles Vermillion, NASA/Calypso Coordinator 
Dr. John Barker, Real Time Landsal Processing 
Mr^ Locke Stuart, Experiment Advisor 
Mr. Albert Whalen, ATS-3 Communications 
* Texas A&M 

Mr. John M. Hill, Oceanographer and Whiting Study Coordinator 
* John Hopkins University 

Mr. Edward Westerfield, Transit and LORAN-C 
Mr. Dan Mitola, Beayondan Coordinator 
* U.S. Coast Guard 

Char les Montanese , LORAN-C 
Robert Riper, LORAN-C 
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^ •Defense Mapping Agency 

Jim Hammack, Historical Data and Navigation Charts 
Dennis Granato, LORAN-C Grid Conversion 

Because of the wide scope of the experiment, separate reports are 
being prepared related to the different aspects of the experiment, 
e.g., real time data processing performed at NASA, Goddard by Dr. 

John Barker and accurate position information obtained by Johns 
Hopkins University under the direction of Mr. Ed Westerfield. This 
report focuses on the remote bathymetry aspects of the experiment. 

A special note of appreciation is given to experiment managers 
Enrico Mercanti and Charles Bohn for their tireless efforts in bring- 
ing together, tmder a tight time schedule, the several elements to ensure 
Che successful completion of the experiment. 

Acknowledgement is given to Charles Vermillion, John M. Hill, 

Rusty Schweickart, and Jack Ford for their support, cooperation and 
participation in the shipboard data collection phase of the experiment. 

Particular thanks to Jacques Cousteau and the crew of the Calypso 
for their dedication and expertise in the underwater phases of the 
experiment, without which, measurements of key experimental parameters 
could not have been made. 

Special thanks to Dan Ml cola and Ed Westerfield and the crew of 
the Beayondan for their patient efforts in obtaining the accurate 
position information, without which, verification of the calculated 
depths with the measured depths would have been difficult. 

At ERIM, the collaboratlcn of Dr. David R. Lyzenga in model 
development and computer analysis of the high gain LandsaC data tapes 
is gratefully acknowledged. 

This report was reprinted after minor corrections were made in 
July 1977. Pages 44 through 127 were repaginated. 
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REMOTE BATHYMETRY USING 
HIGH GAIN LANDSAT DATA 

1 

INTRODUCTION 


This report prepared under contract NAS5-22597 gives the remote 
bathymetry results obtained as part of the NASA/Cousteau Ocean 
Bathymetry Experiment. This experiment using the high gain mode of 
both Landsat 1 and Landsat 2 was conducted between August 21 and 
September 8, 1975 with the test site centered in the Bahama Islands. 
Analysis of Landsat CCT data took place between October 1975 and 
April 1976. 

Two ships, Calypso and Beayondan and thirteen satellites were used 
in the experiment in order to obtain the necessary supporting data and 
position verification. The results are expected to be useful in the 
preparation of charts in order t ) update the location of hazardous 
shoals and to monitor changes in nearshore underwater profiles resulting 
from storm and wave action. 

1.1 SCOPE OF STUDY 

Previously, the International Hydrographic Office defined shoals 
as hazardous to shipping if they fell in the range from 0 to 17 meter 
depths. With the advent of supertankers, the need to "clear*^ the 
shipping lanes to 25 and 30 meter depths becomes a serious concern in 
order to help prevent accidents and thus avoid potential oil spills 
f’"om affecting the environment. More accurate charts could lead to 
savings in ship transportation costs by reducing time spent on unnecessary 
routes to avoid uncertain shoals. 

Not all countries of the world have the necessary ship resources 
to provide the type of measurements needed to update charts in 
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their areas. In order to determine the role of satellite remote 
sensing in helping to provide usable information to alleviate this pro- 
blem, this experiment was conducted to answer two specific objectives: 

1. What is the maximum depth measurable using high gain data from 
Landsat 1 and 2 (in particular, for Band 4 whose wavelengths 
have the best clear water penetration c pability of the 4 
channels available on Landsat 1 and 2)? 

2. What is the accuracy of the depth measurement made using 
satellite data combined with suitable supporting ground data? 

If practical val»’es could be derived from satellite remote bathymetry 
than a comparatively rapid means for surveying the critical ocean areas 
would be created. 

The experiment also included a demonstration of a possible future 
satellite to ship data communication system wherein satellite imagery 
could be relayed directly to a ship n same day of its collection 
giving that ship the most recent information on potential hazards over a 
large area. Since one frame of Landsat data covers a 100 nm square, a 
ship ♦'raveling at 10 knots would be given a maximum of 10 hours of travel 
time over which current information would be available. With two satel- 
lites in orbit, 9 day old information is theoretically achievable for 
any part of the globe, provided a rapid means for processing the data is 
perfected. The real time processing of Landsat data (completed on the 
same day that the image was formed) and transmission of depth contour 
data to the Calypso that was successfully demonstrated during this ex- 
periment is the subject of a separate NASA report. 

If cloud cover lessens the chance for the full realization of this 
concept, the repetitiven'^ss of the satellite coverage does suggest the 
possibility of some type of chart updating at intervals shorter than now 
available. In certain parts of the world charts may be anywhere from 10 
to 100 years old because of the lack of ship and money resources to com- 
plete surveys at shorter intervals. If practical depths are measurable 
at acceptable accuracies, the cost savings in equivalent ship survey time 
could be substantial. 
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1.2 APPROACH 

Remote bathymetry takes advantage of two characteristics. Water 
selectively absorbs different wavelengths of light, and energy at each 
wavelength is strongly absorbed as a function of the depth of the water. 

As the sun’s energy penetrates the ocean, losses occur (1) at the surface, 
(2) through the water column, (3) at the reflection from the bottom, 
then (4) through the water column for the return path, (5) at the 
surface again, and finally (6) through the atmosphere to the satellite 
where it is collected by multispectral scanne* in selected wavelength 
bands. By knowing the sensor sensitivity calibration and taking each 
parameter into account, the MSS voltages expressed digitally can be 
used CO calculate the depth of water. Two variables had to be measured 
at the time of the satellite overflight. The first was the light 
absorption of the water in the same bands in which the satellite sensor 
operates, and the second was the percent reflection of the bottom 
surface. The experimental team aboard the ’’Calypso" made the necessary 
submarine photometer, transmissometer , and reflectance measurements at 
several sites during the satellite overpasses. 

The Calypso made a number of transects to record fathometer depths 
so that satellite-derived depths could be checked against ship measure- 
ments. The team aboard the "Beayondan" operated by the Applied Physics 
Laboratory of John Hopkins University provided accurate location 
information using LORAl^-C and TRANSIT navigation satellite fixes. 

LORAN-C measurements aboard the Calypso were also used for position 
information arJ the ship’s transects were correlated with Beayondan 
measurements at a number of rendevous points. During the time of the 
satellite overpass, both ships were anchored side-by-side so that later, 
the depth measured by the MSS at the corresponding Landsat CCT pixel 
containing the ship could be correlated with actual depth as measured 
by a fathometer. 
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Twelve land points subsequently were used to convert line and 
point elements of the Landsat CCT to geographical coordinates. 

1,3 SUMMARY OF CONCLUSIONS 

Supported only by measurements of average water transmission and 
bottom reflection data, high-gain Landsat Bands 4 and 5 were used to 
construct bathymetry maps of test areas near the Berry Islands and 
near Hollywood, Florida. Depths to 22 m were reliably verified at 
accuracies within 10% (rms) of measured values at the site West of the 
Berry Islands where water transparencies of 0.05 m ^ and 26% bottom 
reflections were encountered. CCT signals two digital counts above 
deep water signals were identified to be caused by light reflected 
from 40 m depths. Landsat data taken in October over the same site 
successfully gave the same depths (within the same accuracies) as the 
depths derived using September data, when both satellite data and 
water characteristics were measured on the same day. 

Depths to 10 m were verified at the Florida site where water was 
less transparent with a = .11 m ^ and bottom reflectance equal to 20%. 

At 3x gain for Landsat 1, Band 4 was found to saturate at 1.3 m 
depths and Band 5 saturated at 0.3 m depth for high reflective bottom 
(about 30% in Band 4). Thus Band 5 can be used to cover the range 
where Band 4 reaches saturation. 

For an area with uniform average characteristics both in time and 
space as seen by the integration over each resolution element, the 
knowledge of point characteristics can be extended to the larger area. 
Alternative ground supporting data such as knowledge of specific water 
depths at two or more control points for each combination of water 
clarity and bottom types would serve equally well. 

Cloud free scenes should be used to avoid spurious shallow 
anomalies resulting from partial cloud cover or from cloud reflections 
of sunlight from the ocean surface. 
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2 

PREVIOUS WORK 

The genera] technique for remote bathymetry was first tested 
using aircraft multispectral data [1]. Two-channel techniques were 
explored to reduce tho dependency on the knowledge of the absolute 
values of bottom refle>t;fon arl water tran^^mission characteristics. 
Two-channel processing was employed in an ERTS-1 experiment [2], In 
that work, the two-channel data gave results to 3 m. At this depth, 
the signal from Band 5, which operates where water more strongly absorbs 
is the limiting factor. Channel 4 data, however, was usable to 9 m 
depths using normal-gain Landsat data. 

The high-gain Landsat data used in this experiment provides a 
better digitization of signal range encountered from sites with 
varying depths. With a wider range of signal spread over more digital 
values, each dig'^tal step represents a finer depth increment so that 
two advantages are achieved. First, depth values can be measured 
within finer bounds with less ambiguity. Second, deeper depths can 
be measured. Two or three digital values above the mean deep water 
background signal (with a standard deviation of less than 2 counts) 
can, in some cases, be interpreted as evidence of the presence of the 
ocean bottom (see Chapter 4). 

Future sensors should be made to operate in bands closer to the 
optimum transmission of water. An investigation of the optimum 
placement for future satellite MSS channels for water depth measure- 
ments is reported in [3]. 
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3 

DESCRIPTION OF EXPERIMENT 

The experimental team met in Nassau, New Providence Island, 
between August 21 and August 25, 1975 to complete prepara* Ions. 

LORAN-C equipment was installed on the Calypso with the assistance of 
Charles Montanese and Robert Riper of the U.S. Coast Guard. Instru- 
ments manufactured by Kahl Scientific Instrument Corporation for 
measuring water transparency and bottom reflections were transferred 
to the Calypso at this time and a trial run took place near the island 
on August 26, 1975 to familiarize the team with the experimental 
procedures. 

3.1 OUTLINE OF EXPERIMENT 

The plan of the experiment called for the two ships. Calypso 
and Beayondan to be on station during the Landsat-2 overpass on four 
consecutive days (August 27, 28, 29, and 30, 1975). The sequence was 
to be reoeated for Landsat-1 on September 5, 6, 7, and 8, 1975. On 
August 27, the first day on station, the rendevous point was north of 
Eleuthera Island. By sailing 90 nm westward during the night (see 
Figure 1) , the two ships were able to be in position the following 
day for a Landsat overpass, west of the Berry Islands on the northern 
edge of the Great Bahama Bank. This site was considered our prime 
test area. It was chosen because of its gradual change in depth from 
one meter of water to deep ocean water in a north and south span of 
25 nm. Laterally to this slope, it is relatively uniform so that a 
broad incline of ocean bottom is available for investigating the 
maximum water depths detectable from a satellite. At the same time 
it enables the depth accuracy to be determined with precision. It 
also makes the exact knowledge of the ship’s position less critical 
during subsequent data analysis efforts, should the position accuracy 
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ORIGINAL PAGE 1S» 
OF POOR QUALTra 


Landsat-2 Aug 30 
Landsat-1 Sept 8 
Center of Pass 4 


Aug 29 
Sept 7 
3 


Aug 28 
Sept 6 
2 


Aug 27 
Sept 5 
1 



FIGURE 1. BAHAMA TEST SITE FOR FOUR CONSECUTIVE LANUSAT PASSES. 
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of a given depth measurement be In doubt within one or two picture elements, 
The water clarity and relatively uniform bottom reflectance at this 
site also permitted the goal of the experiment to be better achieved. 

Table 1 lists the station locations, where underwater measurements 
were taken for each day of the experiment. At each station the 
Calypso and Beayondan were anchored so that accurate location informa- 
tion and depth soundings near each other are well correlated. The 
Beayondan would remain at these stations several hours in order to take 
advantage of the improved position accuracy obtainable from multiple 
passes of the six transit navigation satellites. These satellites 
fixes were then used to improve the LORAN-C position measurements. 

During the time that the Beayondan obtained position information, 
the Calypso sailed at 8 knots along 10 to 15-mile linear transects, 
usually at a compass heading matching the ground track of Landsat. 

Along these transects fathometer data was obtaj.ned continuously with 
minute marks entered on the charts at the time that LORAN-C position 
data were also recorded (see Figure 2) . The transects were also 
chosen so that sometime along the track, the Calypso would pass within 
30 to AO m (one half the size of a Landsat ground resolution element) 
of the anchored Beayondan. The time of passage and the range and 
bearing between the ships were taken into order to maintain correlation 
of depths measured with position (see Appendix A). 

Between stations, the Beayondan also anchored near several 
landpoints to again check satellite position fixes with LORAN-C 
position information obtained from times of arrival differences in 
signals transmitted from a master station and two slave stations. 

Figure 3 shows a typical ship transect defined by LORAN-C. It 
was taken on September 6 starting at 0425 EDT in the prime test site 
west of the Berry Islands. The positions are plotted at one cr two 
minute intervals and were defined by converting the coordinates 
found by the LORAN-C time differences, to geographic coordinates from 
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TABLE 1 

NASA/COUSTEAU OCEAN BATHYMETRY EXPERIMENT 
AUGUST/SEPTEMBER 1975 

TABLE OF STATIOK.S 


DATE 

STATION NUMBER 

LOCATION 

Aug 26 
Aug 27 

A-1 

New Providence Island 
North of Eleuthera Island 

Aug 28 

B-1 

West of Berry Islands 

Aug 28 

B-2 

West of Berry Islands 

Aug 29 

C-1 

Little Issac 

Aug 29 

C-2 

Great Issac 

Aug 29 

G-3 

Great Issac 

Aug 30 

D-1 

Holl)rwood, Florida 

Sept 4 

Transit 

Little Bahama Bank 

Sept 5 

E(5) 

Whiting Study 

Sept 6 

F-1 

West of Berry Islands 

Sept 7 

F-2 

Great Issac 


OF POOR 


FIGURE 2. COUSTEAU AT RECORDING FATHOMETER. Defection corresponds to 
depth while horizontal grid is minute by minute time marks for correlation 
with LORAN-C position information. 



SHIP TRANSECT DEFINED BY LORAN-C 
September 6 Tract Starts 0425 EDT 
and Ends 0517 EDT. 
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maps supplied by the Defense Mapping Agency. The solid line is the 
planned heading. This part of the Atlantic lies at the fringe of the 
LORAN-C network, so that in some cases a poorly locked-on time 
difference gave a displaced apparent position of the ship. However, 
enough reliable points were obtained so that the ships average geographic 
position along the line could be determined. 

3.2 SHIPBOARD OPERATIONS 

A series of underwater measurements was made at six of the 
anchored stations. Three types of instruments were employed to conduct 
these measurements. A multichannel submarine photometer, two filtered 
transmissometers and a combination of two Nikonos underwater cameras 
filtered to match Landsat Bands 4 and 5. The two cameras were used 
to photograph the ocean bottom; a 4~step calibrated reflectance gray 
scale was included in each frame. Three teams carried out the ship- 
board operations to obtain the necessary supporting water data. One 
team operated the multichannel submarine photometer (see Figure 4). 

This instrument consisted of two sensing cells and a deck control 
unit. One sensing cell, called a Deck Cell, remained on board ship 
fastened to the helicopter pad to get an unobstructed view of the sky. 

The Sea Cell, with four identical channels, was deployed away from 
the ship using 200-ft cables and lowered beneath the surface. Readings 
were made from both me.:ers at the surface, and at selected depths of 
1, 5, 10, and 15 m, as well as the bottom. At each level, the divers 
would hold the Sea Cell fixed, it always looked upward, integrating 
the light from the hemisphere above. During this moment, the team 
recording the radiance would cycle through the four detector 
positions. Two positions had filters which matched the Landsat 
Bands 4 and 5; the third was filtered for the blue region; and the 
fourth was unfiltered. 

The necessary range selector adjustments were made each time, 
since less light reached the Sea Cell as it was lowered into deeper 

19 


1 “ 


V 




[ 


I 

i 


I 


formerly WtLLOW RUN LABORATORIES THE UNIVERSITY OF MICHIGAN 



(a) Sea cell being lowered into Zodiac for deployment 200 ft 

away from ship. 






(b) Meters to read light reaching Deck Cell and Sea Cell in 
each of four detector positions. 

FIGURE 4. SUBMARINE MULTICHANNEL PHOTOMETER. 
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water. At the same time that the Sea Cell detector positions were read, 
the corresponding Deck Cell detector readings were also recorded. In 
this way variations in the sun’s illumination during the time between 
readings at different depths could be removed in the data reduction 
phase and only the effect of water depth on the light absorption would 
be determined. 

The photometer data, by measuring the volume- integrated loss of 
energy as a function of water depth, gave values of the water 
extinction coefficient most suitable to the model which was used to 
compute water depth with the satellite CCT data. The. photometer was 
also employed underwater in a second mode (see Figure *>) . By 
restricting its field of view (a tube was attached underwater) the 
diver could obtain readings (again in all four spectral ranges) from 
each step of the calibrated gray scale as well as from typical surrounding 
bottom materials. From these data, the bottom reflectances in Bands 4 
and 5 can be obtained if care in avoiding shadows is taken and no 
bottom sediment is permitted to be in suspension. 

A second technique for accomplishing the same measurement used 
the two filtered Nikonos cameras. A three step gray scale panel with 
measured reflectances of 1%, 10.6%, and 29.6% (see Figure 6) was 
deployed underwater and photographed against the ocean bottom along with 
a seccht disk of known reflectance 84% to provide a fourth calibration 
step. Along with station location codes, the panels contained red and 
green stripes to code the black and white images so that no ambiguity 
would result later in determining in which spectral band an image was 
taken. By measuring the density of the negative across each gray level 
and than at selected points of the ocean bottom, an average percent 
reflectance of the bottom for that particular spectral band could be 
obtained. This parameter was also needed for the calculation of water 
depth. Different mixtures of bottom types were encountered at six 
sites during the experiment. Data from stations F-1 and D-1 were 
subsequently used to make detailed depth measurements. 
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FIGURE 5. DIVER MEASURING LIGHT REFLECTED 
FROM CALIBRATED GRAY SCALE. 


I 

22 



• Mmiti 


2p 

FORMEMLY WILLOW RUN LABOMATORIES THE UNiVtMbiTY OF MICHtOAN 


wuginal page 

OP POOP OT’ • ’ 



FIGURE 6a. THREE STEP CALIBRATED GRAY SCALE MADE 
SPECIALLY FOR EXPERIMENT. Station and Day Codes as 
well as color stripes for coding Nikonos Band 4 and 5 
images were employed. 
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FIGURE 6b. TWO 35 mm LENS NIKONOS CAMERAS FILTERED TO 
MATCH BAND 4 AND BAND 5 OF LANDSAT 1 AND 2. Underwater 
exposures with Tri X film, ASA400 at f/11 at 1/125 sec 
fcr green band and f/5.6 at 1/125 sec for the red band 
were typical. 
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Two transmissometers were also used, each filtered so that beam 
attenuation measurements in spectral bands equivalent to Band 4 and 5 
of Landsat could be made. They were lowered over the side of the 
Calypso (see Figure 7) and the percent light transmission from a light 
source thi ough a fixed distance of water was measured. Measurements 
were made at i, 5, and 10-m depths to check on the assumption of 
uniformity of the water column with respect to the average water 
transmission characteristics. 

Special measurements were made and samples collected at Station E 
in the Little Bah^Mna Bank in order to investigate the source of turbidity 
concentrations that appear in Landsat imagery and are sometimes mis- 
taken for shallow water. The result of this phase of the investigation 
is given in Appendix B. 

3.3 DATA PREPARATION ^ WATER ATTENUATION AND BOTTOM REFLECTANCE 

The results of the submarine photometer work are summarized in 
Figure 8. In this graph, transmission of light, T, where T * e , is 
plotted versus water depth, z. The attenuation coefficient u is usually 
expressed in m Data for both Band 4 (green) and Band 5 (red) are 
given for six stations. The transmission data were de.*-ived from the 
raw photometer data by comparing values in each band near the surfa:e 
with each of the values measured at lower depths. This quotient was 
normalized by referring to the Deck Cell variations recorded at the 
same times as the Sea Cell data were obtained at the sequence of depths. 

The average slopes of each of the lines of Figure 8 were then 
computed to produce the results given in Table 2. From this table, 
tne variations in the type of water encountered at the six sites can 
be seen by an inspection of the values of a for Band 4. The clearest 
water a = 0.052 m occurred at the Eleuthera test site, A-1 and 
theoretically, the deepest depths could be measurable there. The site 
F-1 west of the Berry Islands, measured on September 6, 1975, contained 
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FIGURE 7. TRANSMISSOMETER DEPLOYED Ui^ERWATER TO MEASURE 
BEAM ATTENUATION OF WATER. Distance between light source 
and detector is adjustable and was typically set at 55 to 
70 cm apart. 
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FIGURE 8. GREEN AND RED BAND TRANSMISSION 
VERSUS WATER DEPTH 
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TABLE 2 

AVERAGE WATER ATTENUATION COEFFICIENTS 
DERIVED FROM PHOTOMETER DATA 
BAHAMAS AND FLORIDA 1975 


CALYPSO STATION 

BAND 4 
(m-1) 

BAND 5 
'm-') 

A-1 

Eleuthera 

Aug 27 

0.0522 

0.219 

F-1 

Berry 

Sept 6 

0.0586 

0,314 

B-1 

Berry 

A'jg 28 

0.0638 

0.2"3 

B-2 

Berry 

Aug 28 

0.0661 

0.369 

C-2 

Issac 

Aug 29 

0.0748 

0.326 

D-1 

Florida 

Aug 30 

0.1067 

0.374 
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the next clearest water a = 0.059 m The site D-1 near Florida, 
measured on August 30, 1975, had the poorest water clarity encountered, 
a = 0.107 m 

The experimental test of satellite-measured water depth was 
performed for stations F-1 and D-1 in order to investigate the effect 
of changes in water attenuation on the maximum depth measurable. 
Consequently, the filtered Nikonos imagery of the gray scales for 
sites F“1 and D~1 (see Figure 9) were analyzed to compute the percent 
reflectance of the bottom. The corresponding photometer data taken at 
the same site gave comparable values. 

Transmissometer data in the green band taken at station F-1 
showed that 'he average percent transmission was of 88 + 7 units over 
all depths measured. At D-1, the average percent transmission for 
all depths measured was 78+5 units in the green band. Path was 55 cm 
in both cases. 



FIGURE 9, BAND 4 FILTERED IMAGES. Taken underwater at 
six locations with calibrated gray scale used 
to estimate percent reflectance of bottom. 
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4 

CALCULATION OF WATER DEPTH 


I 


4 • 1 MODEL DEVELOPMENT 

In order to compute water depth using remotely sensed signals, a 
relationship between signal voltage from the Landsat MSS and the depth 
of the water must be developed which includes the parameters that have 
an effect on the sun's energy as it passes through both the air and 
water paths and reflects from the ocean bottom. The complexity of the 
variables in this situation suggest that before a relationship can 
be developed, some simplifying assumptions are necessary. Three 
assumptions were made in this case. First, the model to be used would 
neglect scattering in the water. If we are Investigating maximum 
depth penetration in clear water, this is not unreasonable. Secondly, 
we assume that the sensor signal comes only from direct solar radiation 
defined as E^. For the wavelengths in the green and red bands, this 
again is practical since skylight is dominated by blue wavelengths. 

The third assumption is placed on the p erty of a, the water 
attenuation coefficient, we assume it is independent of the radiance 
distribution for mathematical simplicity. 


In general then, the voltage at the MSS output, V will be related 
to water depth, z, by the equation: 

V - V + V e ■“ ® 

V Vg f Vq e 

**2o(Z 

or ^ ^ small values of 0 and ^ 

where Vg - the signal level from deep water Including atmospheric 

path. 6 and <ti are the refracted solar elevation 
angle and view angle respectively. ^ Is less than 
5® while 6 Is less than 16® 


and 



n 


T 



TT 


31 


I ... 


I 


x: 


” ‘I” 




t 



r 


r 


2pl 


FOMMCRLY WILLOlM RUM L BORATORiCS THE UNIVERSITY OF MICHI<«AN 


such that kg = Landsat-1 MSS 4 sensitivity constant 

^1* ^2 ” surface transmittance (approximately 0.98) 

n - index of refraction = 1.33 

Eq = surface irradiance 

T = e atmospheric radiance transmittance 

and r, = bottom reflectance (assuming Lambertain reflection 
b 

distribution) . 

Hence, if and a are known or can be calculated we can compute 
z by inspecting the Landsat signals after first subtracting the mean 
deep water signal. In a Landsat scene which covers 10,000 sq nm it 
does not prove difficult to find areas of deep water. 

The water attenuation a was measured by the submarine photometer 
and is obtained from Table 2. Most of the terms in are constants 
or known parameters. The bottom reflectance r^, was measured as 
described in section 3. Eq, the surface irradiance was also obtained 
from the Sea Cell reading of the photometer at the surface. E^ could 
also be obtained from knowledge of the sun’s solar irradiance at the 
top of the atmosphere and modified by the atmospheric path transmission 
for a marine atmosphere. Finally, T was assumed to equal 0.8, the 
atmospheric radiance transmittance based on calculations of a Guttman- 
Kwajalein type atmosphere for this site [4]. 


4.2 COMPARISON OF CALCULATED DEPTHS FROM SATELLITE DATA VERSUS 
CALYPSO FATHOMETER READING AT STATION F-1 

On September 6, 1975 at Station F-1 west of the B:jrry Islands, 
in the prime test area, Landsat-1 passed over the Calypso and Beayondan 
at about 10:30 EDT. The computer compatible tapes containing high 
gain data from this scene were subsequently used to construct a depth 
map (see Figure 10). The scene is partly cloudy, but the edge of the 
Great Bahama Bank is clearly delineated, as are several shallower depth 
ranges which were coded by printing in different colors. The digital 




FIGURE 10. BATHYMETRY MAP MADE FROM BAND 4. Site is the northern edge of the Great Bahama Bank 
Clouds or snallow depths less than 4 mare white. Other depths are coded according to legend. 
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values for a given depth range are given in the caption for Figure 10. 
Fortunately, the sky was clear at station F-1 within this scene and 
the satellite sensor signals corresponding to the pixel containing the 
two ships could be identified from the CCT. This was accomplished by 
first rotating and scaling the Landsat data and then developing a 
regression equation that converts line and point numbers from the 
digital format into known geographic coordinates. A twelve-point fit 
using the geographical coordinates of distinct land features in the 
scene were used to compute the regression relations. Since the 
satellite fixes made by the Beayondan gave accurate position information, 
these were used to find the corresponding Landsat line and point 
numbers. The Beayondan’s calculation for the position of station F-1 
was 25° 45.116’ N latitude and 78° 9.001’ W longitude. The corres- 
ponding line and point position in the Landsat data was computed as 
1524, 496 respectively. The digital count for this pixel was found 
to be 66. By careful examination of the digital data it was possible 
to locate which of the six detectors to use in determining the corres- 
ponding deep water signal in an area north of the Great Bahama Bank. 

This value was found to be 59 counts. Using the model described in 
the previous section, the relationship for station F-1 was: 


■ » 


-1 

0.117 


in 


\70.4/ 


where a = 0.058 m ^ was obtained from Table 2. 
Vq = 70.4 for the parameters following: 


Tl = T2 


0.98 


Eq = 12.8 mw cm (Sea Cell reading at the surface) 

T =0.8 (Guttman-Kwajalein atmosphere) 

r^ “ 0.26 (Derived from the underwater photography) 

k = — (NASA-supplied sensitivity constant), 

0.83 mw cm sr 
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The high-gain Band 4 Landsat data for this date saturated at 
127 counts. Substitution into the equation for Z gives: 


^ ■ 0 ^ *” [W] ■ "• 


This is the calculated depth from the assumed model, based on 
careful analysis of the Landsat digital values. 

On board the Calypso, the displacement of the fathometer 
recording pen for this site was 42.1 mm with a scale sensitivity of 
10 ms/20 mm. The temperature of the water at this site was 29.1*C, 
Using the temperature correction of 2.4 m/s/C® above 25®C the velocity 
of sound wao computed to be 1542 m/s. 

The fathometer displacement represents a two-way path for the 
time it takes the sound to reflect from the bottom back to the 
transducer . Thus , 

-3 

o Tc/n-/^ t ^2.1 mm X 10 sec ^ 

Z- = 1542m/s X — = 1542m/s x r = 16.2 m. 

f 2 2 X 20 mm 

The Calypso fathometer transducer is 2.97 m below the water 
line; this bias must be added to the value computed for Z^. Thus, the 
depth Z measured by the Calypso for station F-1 was: 


Z = 16.2 tn +2.97 m = 19.2 m 


The close agreement (2.5% error) between the 19.7 m derived from 
the model and 19.2 m measured by the fathometer gives strong evidence 
to the correctness of the approach and to the potential practical 
utilization of satellite data for water depth measurements. 
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4.3 EXTENSION OF THE CALCULATIONS TO OTHER AREAS 

In order for the satellite technique to be considered practical, 
we must know how well other depths calculated from the satellite MSS 
digital values will compare with ship-measured depths. This necessarily 
raises the question of the reliability of the extension of measured 
values for a and r^ at one site to other nearby areas. This element 
was tested using the ship transect described in Figure 3. 

Along the track, the Calypso position was known from LORAN-C 
measurements and the depth was known for each corresponding minute of 
data samples. A number of positions along this track was selected and 
by similar procedures to that for station F-1 the depths were calculated. 
In addition, the mean deep water signal in Band 4 was calculated for all 
six detectors separately and together. These are listed in Table 3. 

One can see that the Landsat detectors are not equally noisy so 
that depth accuracy can vary from point to point in a scene. For this 
analysis, the standard deviation of 1.77 calculated for all six 
detectors was chosen to define the error bars on the depth calculations 
from the satellite data. A graph showing the calculated values versus 
the ship measured depths is given in Figure 11. 

The ideal case would be for all values to fall on the straight 
line. It can be seen that agreement to 22 m depths is readily 
acceptable but beyond this range the calculated values are consistently 
lower than measured values. It is encouraging to note that for this 
site, depths at 40 m were detectable (even though they were not 
correctly measured as such) . There could be several reasons for the 
short estimate. 

First, the model has neglected scattering effects. At longer 
water path lengths, this assumption could be invalid. Second at deeper 
depths, less bottom vegetation may occur, so that the effective 
reflection could be higher and thus appear to be a signal from a 
shallower depth. In any event, the signal from the 40 m depth represents 
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TABLE 3 


AVERAGE MSS4 SIGNALS OBSERVED 
OVER DEEP WATER, 

SEPT. 6, 1975 


Detector 

Mean Signal 

Std. Deviation 

1 

57.80 

1.43 

2 

57.34 

2.05 

3 

57.85 

1.99 

4 

58.68 

1.55 

5 

57.70 

1.79 

6 

59.09 

1.05 

all 

58.14 

1.77 


Based on sample of 2100 Points for 
each detector. 
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(Calypso fathometer - meters) 


FIGURE 11. 

LANDSAT-1 CALCULATED DEPTH 
vs MEASURED DEPTH 
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only two counts above the deep water signal, indicating that the 
effect may also be statistical. Note the location of the upper eiror 
bar on the 40 m test case. It is close to 35 m. An error of 5 ro out 
of 40 m is 13%. The agreement is better for shallower depths where 
the signal- to-noise count is more favorable. 

4.4 EXTENSION OF THE CALCULATIONS TO ANOTHER DATE 

A stronger test of the metliod can be made by assuming that the 
ground data taken in September is stationary in time and using it to 
calculate depths with Landsat CCT’s obtained in October. On 
October 12, 1975, Landsat-1 obtained coverage of the station F-1 area. 
Since the solar zenith angle was different, the signal levels were 
found to be lower. However, the noise variation as seen in the 
calculation for the standard deviation in Table 4 were similar to the 
September case. For October, the mean signal had dropped to 51.1 
counts and the standard deviation was 1.63 counts. The signal change 
was proportional to the cos 0 where G is the solar zenith angle. 

In a similar fashion as above, a graph comparing Landsat cal- 
culated depths versus measured depths was constructed (see Figure 12) . 

The water transparency value and the bottom reflection value used in 
the model were those measured in September. 

The October satellite-calculated depths followed a similar pattern. 
Good agreement in the 20 to 23 m range with the error on the short 
side for the deeper depths. lnter-«9tingly , the upper ?rror bar for 
this date for the 40 m data point now cresses the correlation line. 

If the signal-to-noise of the Landsat signal could be improved 
even by two counts better reliability woula be achieved for calculated 
measurements to 40 m. 

The results shown in Figures 11 and 12 suggest that remote 
bathymetry from space is nracical. If the measurements of water 
clarity and bottom reflectance at a small number of control points can 
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TABLE 4 

AVERAGE MSS-4 SIGNALS OVER DEEP WATER 
OCT 12, 1975 LANDSAT -1 

Detector Mean Signal Std. Deviation 


1 

51.11 

1.21 

2 

51.45 

1.57 

3 

51.08 

1.49 

4 

50.52 

1.93 

5 

51.19 

1.36 

6 

51.60 

1.89 

all 

51.15 

1.63 


Based on sample of 2100 points for each detector. 





MEASURED DEPTH (meters) 
(from Sept. 6, 1975 fathometer data) 


FIGURE 12 

LANDSAT-1 CALCULATED DEPTH 
BASED ON OCTOBER 12, 1975 DATA 
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be made and if they represent values that can be used over a few months 
in time, then updated charts for the locaLion of water depths hazar^’^us 
to ships is feasible. And not only can position informatior be derived 
from space data, but estimated depths with accuracies on the order of 
10% or less in the range from 0 to 22 m is feasible with present 
satellites in orbit operating in the high gain mode. Suitable supporting 
cata is necessary for highest accuracy and can take several forms. 
Knowledge of the depths at selected points can serve the same role as 
knowledge of water uiari y and bottom reflectance. However, there is 
a world wide distribution of measurements from different countries of 
water transparency and bottom reflectance data that could be obtained 
and these values tested for determination of overall accuracy. One 
advantage of the 70 m resolution of the Landsat sensor is that small 
variations in the bottom surface are averaged out. With a library of 
parameters for different areas and with careful analysis of repetitive 
satellite data, charts of ocean areas with potential hazards can be 
constructed. 

To illustrate what is possible in map construction, a relatively 
cloud free area of the October 12, 1975 scene for the area west of the 
Berry Islands was used to make a depth map. The computer output of the 
CCT was first rotated and scaled in order to produce a final product 
with proper aspect at 1:300,000, The scale was chosen to match chart 
N.O. 26320 covering the Northwest and Northeast Providence channels. 

The area on this map of the investigation is labeled Northwest channel 
with a reference to a very intricate passage due to hazardous shoals. 

The first edition of the chart is given as June 7, 1965 and it was 
was revised February 12, 1973. A transparency of this mao f'^r the 
area selected from the Landsat CCT was used as an overlay to the color 
coded depth map generated by ERIM using the color iniv-jet printer on 
the output of the MIDAS computer, MIDAS is a multivariate interactive 
digital analysis system in developmental stages under NASA sponsorship. 
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The map made from Band 4 is illustrated in Figure 13. The color 
association for each range of depth is given in the legend. The numbers 
on the overlay express the depth in fathoms (6 ft or 1.83 m) . 

The effect of noise fluctuations in the data can be seen by the 
texture of the dark blue area representing dee:- water. This is relaced 
to the differences in the mean signal levels and their .ariances for the 
six detectors as shown in Table 4. 

The spurious shallow depths seen in deep water areas in the lower 
right are the result of the presence of clouds, cloud shadows and sun- 
light reflection from clouds onto the oceaii surface. Obviously these 
artifacts can be eliminated by choosing only cloud free scenes or 
cloud free areas within a given scene. The ge..jral agreement between 
the overlay and the satellite map is encouraging. For example the pink 
yellow boundary which represents about 5 m depths falls on the chart 
between 2 to 2.5 fathoms or near 3.7 to 4.5 m. 

One method to reduce the effect of the difference in signals from 
the six detectors operating in one band is to average their values line 
by line by a six by six **moving window” algorithm operating on the CCT 
data. This procedure was implemented and the results are shown in 
Figure 14. The procedure smooths the data so that each color’s depth 
range is less ambiguous. The trade off in producing this type of map is 
that the spatial resolution is made coarser. The position uncertaincy 
has been increased. However, position information can be obtained 
before the smoothing process. 

Another benefit for depth charting is that with a 6 x 6 element 
smoothed data set, a single digital value can be printed which will 
represent a narrow depth range and produce a contour line depth profile 
map as shown in Figure 15. 

Each ^ ine represents the depth boundary between adjacent depth 
ranges of Figure 14. Where the line widens that area would be inter- 
preted as having the same water depth. Other forms of CCT processing 
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FIGURE 13. BATHYMETRY MAP FROM LANDSAT-1, BAND 4, OCTOBER 12, 1975. 
With overlay from Chart N.O. 26320. 
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FIGURE 14. BATHYMETRY MAP FROM LANDSAT-1, BAND 4, OCTOBER 12, 1975. 
With overlay from Chart N.O. 26230; Data averaged over 6X6 resolution 

elements . 
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are possible using the satellite data and further research could produce 
those output formats most helpful to the map using community. 

4.5 SATURATION OF BANDS 4 AND 5 

As seen from Figures 13 and 14, the red arer‘=- represent depths not 
measured because of signal saturation. With the nigh gain mode for Landsat, 
there will be an increase in the areas where signal saturation will occur 
in Band 4. This occurs in shallow depth areas, with high reflective 
bottom materials on the order of 30% or higher. 

This aspect was further investigated for an area just north of 
Andros Island. Bright calcareous sands almost awash produce signal 
levels at 127 counts. A plot of the signal level difference versus water 
depth is given in Figure 16 for both Landsat Bands 4 and 5. The signal 
level associated with deep water for each band was subtracted from the 
Landsat signal to give the linear relation shown on a semilog graph. 

For this case depths were estimated from the chart but the error bars 
are calculated from the deep water noise fluctuations as before. 

The data shows that for this area with the conditions prevailing, 

Band 4 saturates at 1.3 m while Band 5 saturates at 0.3 m. Figure 16 
also explains why two channel processing for mapping depths with 
present Landsat MSS is limited. The minimum and maximum depths where 
both channels have usable data are from between 1.3 m and 6 m. This 
could be improved in the future by two channels operating in the blue- 
green and green bands of the spectrum. 

The use of Band 5 (also 6 and 7, if needed) to isolate water depth 
at small intervals can bo seen from Figure 17. There the amount of area 
(red color) giving a saturated signal is greatly reduced. The finer 
depth increments show the complimentary nature of Band 5 with regard to 
the range where Band 4 saturates. The higher absorption of radiation in 
Band 5 helps to Improve the depth division with steps fraction of meters 
apart. In surface water hydrology or shall ow-water barge-transport 
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WATER DEPTH (meters) 


FIGURE 16. SATURATION DEPTHS FOR LANDSAT-1 BANDS A AND 5 


51 



1 

] 



— A. 










{ 

‘ff 

i 

1 


problems, such narrow level mapping of depths should prove 
advantageous. In a similar manner. Band 6 could be used to "slice" 
the shallow depth ranges where Band 5 saturates. 

One notable imprr vement in shallow water mapping with Band 5 is 
that a clear passage of 3.6 m depth is shown for the Northwest channel 
(lower middle of Figure 17). In Figure 13, for the same area, a red 
zone defining a range of 0 to 5 m leaves a clear passage ambiguous. 
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5 

DEPTH CALCULATIONS NEAR THE FLORIDA COASTLINE 

We have seen the results for remote bathymetry from satellite data 
for a test area with clear water characteristics, that is, when the water 
extinction coefficient was equal to 0.058 m On the fourth day of 
the experiment (August 30, 1975) uhe Calypso anchored near the coast of 
Florida and obtained data at station D-1, to determine the loss of light 
penetration in areas of different clarity (see Table 2) , 

The satellite CCT data for this date was not available but the 
site was imaged on September 7, 1975 under partially cloudy conditions. 

A portion of the east co^^st of Florida, extending approximately 
from Lake Worth to Biscayne Bay, is contained within Landsat frame 
5141-14503 (September 7, 1975). A subset of this data, covering 
the shoreline from Port Everglades to Miami Beach, was rotated and scaled 
in order to yield a geometrically correct display on the MIDAS ink-jet 
printer at a sc'le of 1:80,000. The scale of the transformed data set 
was verified by means of an analysis of four ground control points, and 
the following relationships were established for the line and point 
numbers on the rotated data: 

line no. = 1871 - 28.06 x lat 
line point= 280 - 23.88 x long 

where lat is the latitude in minutes north of 25°N, and long is the 
longitude in minutes west of 80®W. 

The Beayondan position information from the satellite inter- 
rogations gave the coordinates for D-1 at 26® 3.8556’N and 80° 5.648’W. 
Using the above relationships, the line and point numbers for Station 
D-1 were determined to be 79 and 145, respectively. A display of the 
data from the CCT showed the signal values at this location to be 68 
counts in MSS 4 and 34 counts in MSS 5. 
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An area of Jeep water further offshore was found to have a mean 
signal of 61.8? vO’.nt;-: with a standard deviation of 2.99 counts in 
MSS 4, and a mean signal in MSS 5 of 34.61 counts with a standard 
deviation of 2.84 counts. 

The signal variation in the deep water area for this data was 
found to be higher than that found near the Berry Islands. We inter- 
pret this to mean that the wat-^r color (and hence the returned signal to 
the sensor) near Florida is considerably more variable spatially 
because of the presence of the Gulf Stream, river discharges and ocean 
outfalls, and other industrial effluents along the Florida coastline 
which contribute to the variation in light scattering and absorption. 

The signal level of 68 counts for station D-1 against a mean deep 
water signal of 61.8 counts clearly indicates that a bottom reflected 
component exists in MSS Band 4. This is not true for Band 5, where 34 
counts at station D-1 is less then 34.6 counts for the deep water 
signal. 

In order to calculate the depth at station D-1 the various ground 

measurements, Including the underwater photography of the gray scale 

against the bottom, were analyzed to yield a bottom reflectance of 

20%; the photometer data gave an extinction coefficient of 0.1067 m ^ 

in MSS 4 at this location. As a practical matter, the film negatives 

with the gray scale are scanned with a densitometer and the film 

density (or traismission) is plotted versus the known panel reflectance. 

An average film density (or transmission) was measured for the bottom 

and the bottom reflectance is obtained. However, this quantity 

measureu from data taken under water can be defined as the wet 
^b 

reflectance = where was defined as bottom reflectance in 
Section 4 and n is the index of refraction. When we report r^ = 20% 
for D-1 or r^^ = 26% for F-1 these values are the dry reflectances. 
Substituting the appropriate values in the depth equation (see Section 
4) gives the following relationship between water depth and MSS 4 
s ignal : 
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= 54 for this site since the r. is different from that 
^ ° of site F-1 


V - 62 + 54 e 
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This equation yields i depth of 10.2 m at Station D-1. 

The fathometer reading for this site on August 30 gave a depth 
value of 11.2 :r. following the procedure described for site F-1. Again 
the agreement is encouraging. The higher extinction coefficient, 
however, limits the maximum depth measurable at site D-1. In order 
to determine the range of depth detectable along the Floriaa shore- 
line, two transects perpendicular to shore were investigated. One 
transect covered the track of the Calypso for this site and is shown i.. 
Figure 18a. The second ..ransect was taken near Sunny Isles (25° 59.75’N) 
and corresponded to a fathometer line measured two years earlier by 
ERIM personnel in the same test area. This is shown in Figure 18b, 

The plot shows CCT count versus point number or position offshore. The 
CCT data was averaged in rectangular arrays 5x1 pixels with the 1 ^ng 
axis parallel to the cof 'tline. The fairly smooth curves can be used 
to investigate the maximum depth penetration under these more difficult 
circumstances. Figure 19a and b show the conversion to depths using the 
equation given above. 

For the Sunny Isles transect, a plot of the ground truth is also 
shown for comparison. Good agreement is found to depths of 5 m within 
the first 1,000 m from ohore. Subsequenc points further from shore 
were calculated to be deeper than measured depths. This implies a 
change (lower) in the value for the bottom reflecMon. A corrected value 
would have to be Inserted into the model for this part of the transect. 
Beyond 2,200 m from shore, where actual depths are greater than 14 m, 
the satellite depths are shallower than measured values, producing the 
effect noted for clear water beyond a depth of 22 m (see Chapter 4). 
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A color coded depth chart for this site was aJso produced on the MID\S 
Ink-jet printer (see Figure 20). Depths from 0 to 15 m are indicated 
in 3 m increments by means of six colors op this chart. These colors 
and the corresponding data and depth ranges are indicated in Table 5. 
Land is indicated by overprinting in black those areas where the 
signal in MSS 7 is greater than 12 counts. A few pixels are also 
visible with a pink color: these indicate pixels where the MSS 4 

signal was saturated in high gain mode. The data presented in this 
display were smoothed 5 lines by 1 point in order to reduce random 
noise while maintaining spatial resolution in the perpendicular-to- 
shore direction. 

A cloud present in the image nearSiiore gives the anomalous 
shallow water reading. The scale of 1:80,000 was matched to chart 
C&GS 11466 from which the overlay to the color coded map was made. 

Numbers of the chart are expressed in feet. ihe correspcnden^e 
of the location of underwater variations along shore ia the chart to 
that of the satellite map is favorable. The feasibility of the 
technique is again demonstrated even though only shallow depths are 
measurable because of the higher extinction coefficient. The maximum 
depth discernable along the track containing station D-1 is about 18.5m. 
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TABLE 5 

COLOR CODE FOR DEPTH RANGE 
OF FLORIDA COASTLINE MAP 
DEPTH 
DATA RANGE 
RANGE (meters) COLOR 


1- 

64 

> 15 

Purple 

65- 

66 

12-15 

Light 

Blue 

67- 

70 

9-12 

Green 

71- 

77 

6- 9 

Tan 

78- 

90 

3- 6 

Yellow 

91- 

126 

0- 3 

Red 

Saturation 


Pink 



Land 

Black 


FIGURE 20. 

COLOR CODED DEPTH MAP 
MADE FROM LANDSAT IMAGE. 
Taken September 7, 1975. 
Overlay Is from Chart 
C&GS 11466 at approxi- 
mately 1:80,000; sound- 
ings are In feet. The 
anomalies offshore are 
due to clouds. 
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CONCTAISTONS 

The success of the satellite calculated water depth using high 
gain Lcmdsat MSS 4 data demonstrated the feasibility of remote 
bathymetry from space. In clear water (extinction coefficient = 0.058m 
depths to 22 m can be r.easured reliabl}^ and signals from 40 m depths 
can be differentiated sometimes from deep water signals that are due to 
scattering only. In less clear water (extinction coefficient = .11 m 
depths to 10 m are reliably measured and detection of the presence of 
18,5 m depth is feasible. For these values to be generally realizable, 
the bottom reflection should characteristically be between 20 to 26% in 
the green band (Landsat MSS 4), The best accuracy of the depth measure- 
ment will occur if knowledge of the areas average water transmission 
characteristics are known as well as the reflectance of the bottom. 

Single-cha*:nel charting is reliable where assumption of uniform 
bottom reflectijn and water transmission characteristics are valid. If 
changes do occur, adjustments of the appropriate model input parameters 
t mur ' be nade to maintain accuracy. Otherwise those pixels should be 

edited. The technique may use alternative information to bottom 
! reflectance and water clarity. Knowledge of depths at control points 

} within thw scene, made available by water level gauges or by aircraft 

I with laser depth- ranging equipment and ship support measurements are 

I equivalent, A catalog of average values for an area will give useful 

! a p^“i( ri information. Repetitive satellite coverage helps provide 

‘ infovination as to possible deviations from average values. 

I Under certain conditions of bright bottom reflectances, 30% or 

I . greater and shallow depths, less than 1,3 m, the high gain data of 

'j 

j Band 4 will saturate. In this case, it is possible to use Band 5 data 

I to determine depth in fractions of meter steps for the range where Band 4 

I saturates because of the greater light attenuation in Band 5. 




2pi 


FORMEMLY WILLOW RUN LABORATORIES THE UNIVERSITY OF MICHIGAN 


Cloud free Landsat scenes offer the best solution to depth map 
construction to avoid anomalous shallow depths due to clouds and ocean 
reflection of sunlit clouds. Incorporation of scattering effects into 
model relating Landsat signals to water depth could improve accuracies 
of deep water estimates. 

Detection of deeper depths or improved reliability of intermediate 
depth calculations would be possible if optimum band location could be 
employed on future spaceborne sensors. An analysis of the effect on 
S/N ratio by placement of the channel bandwidth between .45 to .60 pm 
was made on 0.1 ym intervals [3]. That study concluded that improve- 
ment in S/N would occur by placement of the band in the range 
.47 to .57 pm. The same methodology of that study was employed to 
determine placement with .07 pm bandwidth intervals. The results 
are shown in Figure 21. The analysis was done for three types of 
water clarity and two atmospheric conditiou.^. The graph shows the 
improvement to be gained by proper location of uit channel bandwidths. 
For clear water, 0.45 to 0.52 pm would be optimum. As a compromise 
to all three types of water, 0.47 to 0.54 pm might be selected. 

The NASA/Cousteau Ocean Bathymetry experiment also demonstrated the 
potential of near-real-time transmission of processed satellite data 
to ocean going vessels through ATS-3 communication link. The report 
of this phase of the experiment is given in Reference 5. 
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APPENDIX A 

Technical Memorandum S3R-76-067 

ANALYSIS OF LORAN-C DATA COLLECTED ABOARD 
CALYPSO DURING OCEAN BATHYMETRY EXPERIMENT 

April 1976 

The Johns Hopkins University 
Applied Physics Laboratory 



the JOHNS HOTKINS UNIVERSITY 

APPLIED PHYSICS LABORATORY 

LAund marviano 


S3R-76-067 
IW10S3R0 
April 19, 1976 


TO: Distribution 

FROM: E. E. Westerfield 

SUBJECT: Analysis of LORAN-C Data Collected Aboard Calypso 

During Ocean Bathymetry Experiment 


BACKGROUND 

In July of 1975 NASA requested that APL investigate 
the feasibility of accurately determining the position of Jacques 
Cousteau's ship Calypso during a joint Cousteau Society/NASA 
operation. The purpose of the operation was to determine the 
feasibility of accurately determining depths of water in shallow 
areas utilizing data from the LANDSAT satellite. APL was speci- 
fically requested to look into installing a TRANSIT Navigation 
receiver aboard the Calypso and to work with the Coast Guard in 
the use of LOPJ\N-C. Preliminary studies indicated that it would 
be difficult to install TRANSIT gear aboard the Calypso due to 
very limited space, plus an uncertain power system. This was par- 
ticularly true since Backpack equipment, which is capable of 
battA^ry operation, was not in any condition to be utilized at that 
time and time was not available to put it in condition. In ad- 
dition, it appeared that Calypso would be constantly in motion 
and that it would be difficult to obtain a sufficient number of 
passes at any one point for high quality position determination. 

It was suggested by Laboratory personnel that the TRANSIT 
gear be installed aboard another vessel that could support l 
gear adequately and that would accompany the Calypso. The Labi *a- 
tory's 85 ft. Motor Sailboard Beayondan was picked for this job, 
principally because it was available and required a relatively 
small crew to operate. 


**rj'x;e5ding page blanit not fiuod 


I 


r 


THE X)HN8 HOPKINS UNIVERSlTV 

APPLIED PHYSICS LABORATORY 

LAURIi. Maryland 


S3R-76-067 
Page -2- 


The proposed method of operation was for the Beayondan 
to sail to one or two sites in each area and to anchor long 
enough to accurately survey the location by means of the Navy 
Navigation Satellite system (TRANSIT) . The Beayondan was also 
equipped with LORAN-C gear so that the LORAN grid could be ac- 
curately calibrated at each point. 

At NASA's request the Coast Guard installed a LORAN-C 
receiver aboard the Calypso. The receiver was a unit that APL 
had purchased for the Coast Guard and loaned back to the Coast 
Guard for this operation. Coast Guard personnel installed it 
aboard the Calypso In Nassau. This receiver was a Deca ADL81. 

The purpose of this receiver was to determine Calypso's position 
while it was underway and too far from Beayondan for Calypso's 
position relative to Beayondan to be accuractely determined 
using visual means. 

PURPOSE OF THIS REPORT 

The purpose of this report is to serve as a vehicle to 
disseminate the data recovered by processing the LORAN-C data 
collected aboard Calypso, during four of the data collecting runs 

OVERVIEW OF LORAN SYSTEM 


LORAN is a hyperbolic radio navigation system that makes 
use of a master and two or more slave stations Each station 
transmits a short pulse of lOOKC carrier at a rate of approximately 
10 times per second. The exact rate varies from chain to chain. 

All LORAN-C stations transmit on a frequency of 100 KHz. Precise 
navigation cannot be obtained by using the pulses, so a system is 
utilized whereby measurements are made using the 100 KHz carrier. 
This is done by a method known as cycle matching. The third posi- 
tive going zero crossing of each pulse is considered the epic 
point for the pulse. 
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The LORAN-C receiver is built in such a manner that 
it can track a master station and up to three slave stations 
at once. The receiver displays the difference in arrival time 
of the signal received from the master station and each of 
the (TD) slave stations being tracked in microseconds. The re- 
ceiver used aboard Calypso was capable of tracking three slaves. 
One of the difficulties with the receivers and the system in 
general, is that it is necessary that the receiver lock on and 
track the third cycle as received from tb'^ master and all the 
slaves. When the signals are very weak the receivers tend to 
track the wrong cycles, resulting in an error in the time dif- 
ference measurement, i.e., the difference in time between when 
the master station was received and a particular slave, by some 
multiple of 10 microseconds. This can result in a significant 
error in navigation. 

In the area of operation the master station, which io 
located at Cape Fear, is received quite well, as well as the 
signal from Jupiter, one of the four slave stations of this 
chain. The signals from the slaves located at Dana, Indiana, 
Newfoundland and Nantucket, however are quite weak. This re- 
sulted in the receiver often tracking on the wrong cycle or 
switching cycles in the middle of the data run. In addition, 
the geometry of the LORAN-C stations as viewed from a user in 
the area is poor. The lines of position (lines of constant time 
difference) obtained between the master and the various slaves 
are running very near parallel. As a result, a very small error 
in measuring the time difference can result in a relatively 
large position error. This is particularly true as one approaches 
the coast of Florida, where one is near the baseline of the 
Cape Fear - Jupiter pair. 
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In addition, It appeared that the LORAN-C receiver 
aboard the Calypso was not functioning too well. The receivers 
are equipped with a series of lights to indicate when they are 
properly locked on the stations. This receiver often gave an 
indication of being properly locked even though it was not locked 
or was locked onto the wrong cycle. 

DATA COLLECTION 

The data was manually recorded by Calypso personnel 
on log sheets during the various data collection runs of Calypso. 

The runs were arranged so that typically Calypso ran very near, 
i.e., usually within 50 feet of Beayondan, at least once during 
the run. This proved to be invaluable in that the position of 
Beayondan, which had been obtained by means of the Navy naviga- 
tion satellite, could be utilized to calibrate the Calypso LORAN-C 
receiver. This information could be used to determine which 
stations were being properly tracked, i.e. , if the receiver was 
tracking on the correct cycle and if not what the correction factor 
was . 

PROCESSING OF DATA 

To obtain the best estimate of Calypso position while 
it was underway, was a very involved process due to the relative 
poor quality of the data, pxus occasional errors in the hand re- 
cording process. 

The data was first punched on IBM cards. The cards were 
then read into a Honeywell H21 computer where they were stored on 
disc. Use was then made of a Hewlett-Packard 9821 calculator 
with associated plotter. Hardware was available that allowed the 
data to be read directly from the disc and entered into the cal- 
culator. The first step in the processing was to plot the time 
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difference that haa been read from the card against time. Jumps 
of 10 microseconds in the data could be readily seen from the 
plots plus other errors such as punch errors could be readily 
seen. The cards were corrected ui'til a smooth curve was ob- 
tained from each set of TD's. 

The TD’s collected by the Calypso receiver at the time 
of closest approach were compared to TD’s computed for the posi- 
tion. Routines were then utilized for offsetting the TD’s 
normally by a multiple of 10 microseconds, to make them nearly 
agree with the computed TD’s for the point of closest approach. 

The TD’s for the complete run were then smoothed by 
utilizing a linear regression algorithm where a straight line is 
fitted to each set of 7 consecutive points, in a least squares 
process. The equation of the line is then used to calculate the 
TD of the center of the line. This was repeated for each data 
point, always utilizing the data collected for that point and 
three points on each side. In the case of the end points, the 
value is calculated from the coefficients of the nearest 7 point 
segments . 

Utilizing the smooth TD’s, the position at each point 
was calculated utili ing a standard LORAN-C computational algorithm. 
The positions were then plotted along with the positions of both 
Beayondan and Calypso when at anchor. In some cases the data was 
still quite ragged and additional smoothing was done directly on 
the latitude and longitude number utilizing the same algorithms 
used for the TD’s, i.e., where a straight line is fitted to 7 
consectuive data points. 

prcx:essing of data 

Four sets of data collected aboard Calypso have been 
received from Mr. Fabian Polcyn of the Environmental Research 
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Institute of Michigan, who was principal investigator for 
the experiment. This data covers the following time spans. 


27 

August 

04:30 

to 

06:48 

local 

time 

27 

August 

10:36 

to 

11:23 

local 

time 

28 

August 

03:17 

to 

05:48 

local 

time 

6 

Sept. 

04:25 

to 

07:57 

local 

time 


The unique features of processing each data run is discussed 
below. 

27 August 04:30 through 06:48 

During the majority of this run personnel aboard 
Calypso recorded data from three slave stations. All three 
possible combinations were then processed, i.e., Jupiter/ 
Nantucket; Jupiter/Newfouudland and Newfoundland/Nantucket. 
Following the position computation, data from each station pair 
was plotted. As was predictable knowing the typical strength 
of the signal, the Jupiter/Newfoundland data was by far the 
best, being very smooth and not requiring any post position 
computation filtering. The other data was much noisier, with 
the Newfoundland/Nantucket data being essentially useless. 

This was due not only to the fact that the signals were weak, 
but also that the paths to both slave stations from the vessel 
are nearly coincident and as a result a very poor geometric 
solution is obtained. Data from this pair deviated by over 4 
miles at times from the data of the other pair. The difference 
between Jupiter/Newfoundland and Jupiter/Nantucket was half a 
mile at the worse. The Jupiter/Newfoundland data agreed much 
better with the NAVSAT data. Figure 1 is a plot showing the 
results of computations. utilizing all thvee pair. Figure 2 
is a more detailed plot showing only the Jupiter/Newfoundland 
data. X's are shown at each data point. Data was typically 
collected every minute and the time is ini icated every 10th point. 
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In addition, the positions of the Beayondan at six o'clock and 
at eight o'clock are plotted as well as the position of the 
Calypso before and after the data runs. The position of the 
Beayondan was obtained via the Navy navigation system and is 
accurate to around 15 meters for the 8 o'clock position and 
to 100 meters for the 6 o'clock position. The error of this 
fix is due to the fact that only one pass as taken so that no 
exact statistics can be arrived at. Figure 2 not only shows 
the current run under discussion, but also the run made between 
10:26 and 11:23 local. Figuj^e 3 shows the UA.ta for this run 
without the time marks. The turn at the end of the run can be 
seen from this plot. Table 1 is a listing of the processed 
data giving the best estimate of the latitude and longitude at 
each time point. For this run Jupiter and Newfoundland re- 
quired essentially no correction at all. The run made \itiliz- 
ing Newfoundland required a 10 microsecond correction, i.«;., 
the receiver was 1 cycle away from the proper third cycle or 
tracking. 

27 August 10:34 to 11:30 

Data for this run was available from three different 
slaves. As usual the data from the Jupiter slave was good while 
that from Nantucket and Newfoundland was weaker with the New- 
foundland/Nantucket pair oeing weak enough that it .-.s not 
utilized. In the case of Newfoundland the receiver was 4 cycles 
off, i.e., 40 microseconds. This error was detected by compar- 
ing the Calypso TD’s to that collected aboard the Beayondan. 

The actual comparison was done when Calypso was anchored and 
range and bearing was taken from Beayondan to Calypso via radar 
and hand compass. The Beayondan position was from NAVSAT. 

Figure 4 is a chart showing the data and Table 2 lists the pro- 
cessed data. The position of Beayondan and Calypso when at anchor 
is again shown on the chart. The reader should note the high 
scale factor for this particular chart. 
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28 Augus^ 03:19 to 05:32 

Data was again recor'led from Jupiter, Nantucket and 
Newfour.dland. Again the signal from Nantucket was obviously 
very weak as the receiver was often unlocked. All data pro- 
cessing was therefore done with the signa?r from Jupiter and 
Newfoundland. The signal however, was off by 60 microseconds, 
i.e., 6 cycles, which was corrected in processing. The time 
differences were filtered prior to the position computation. 
Figure 5 is a plot of the computer positions prior to the in- 
sertion of post processing filtering, while Figure 6 is a 
plot of the position after filtering. The position of the 
Beayondan at its anchor point is also plotted. It appears that 
the second filtering probably improves the data becfuse it is 
very unlikely that the ships made sudden changes of course as 
indicated by the unfiltered plot. Table 3 is a listing of the 
filtered positions vs time. 

6 September 04:25 to 07:57 

Data was collected from the three slave stations. 
However, for this run the Newfoundland data was very sporadic, 
while the Nantucke .. data was quite usable. Nantucket/ Jupiter 
was therefore used for all data runs. No major corrections were 
required, i.e., the receiver appeared to be tracking on the 
proper cycle. The data was sufficiently noisy, however, that 
post processing filtering was utilized. Figure 7 is a chart 
showing movement of Calypso while Table 4 lists the position as 
a function of time. Tables 5 and 6 summarize the results of 
'TWSAT data collected aboard Beayondan, including the bias errors 
noted for the LORAN-C grid. The position computed for Calypso 
when at anchor is also shown. This was computed using Beayondan 
position and offsetting using range from either radar or optical 
measurement devices, and bee ring obtained by hand compass. 
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CONCLUSIO N 

The LORAN-C results shown here are probably accurate 
to around 200 meters, i.e., 1/10 of a mile. This number cannot 
be proven but it is based on the accuracy of closures with 
Beayondan typically found in processing the data. Considering 
the difficulty with the LORAN data.- the availability of Navy 
navigation satellite data made it possible to process the LORAN 
data with reasonable accuracy. LORAN-C and TRANSIT are an 
excellent combination, in that the TRANSIT can give high accuracy 
when the ship is at anchor while LORAN-C can provide reasonable 
accuracy when underway. It is felt that if the operatic*! had 
been carried out in an area where LORAN-C was better, then the 
data processing operation would have been simpler and the re- 
sults would have been obtained much faster. 
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Table 1 




TIME 


lat; ri)D£ 

longitude 

U30 

?5 

3T. 99900 

•76 

52.63158 

u31 

25 

33.9228(» 

-76 

52.49940 

432 

25 

33 ,841 98 

-76 

52.36596 

433 

23 

33.76025 

-76 

52.24332 

434 

25 

33.69024 

-76 

52.12248 

43S 

25 

33.59640 

-76 

51 .99294 

436 

25 

33.50580 

-76 

51 .87030 

437 

25 

33.40428 

-76 

51 .75090 

438 

?5 

33.32604 

-76 

51 .65346 

439 

25 

33.22800 

-76 

51 .54348 

440 

25 

33.1 4034 

-76 

51 .44358 

441 

25 

33.05868 

-76 

51 .33672 

442 

25 

3.?. 99706 

-76 

51 .22652 

443 

25 

32.96436 

-76 

51 .1 4736 

444 

25 

32.94324 

-76 

51 .06990 

445 

25 

32.95290 

-76 

51.001 80 

44o 

25 

33.001 38 

-76 

50.95374 

U 47 

25 

33.05988 

-76 

50.91 324 

U48 

25 

33.1 4532 

-76 

50.88840 

449 

25 

33.25884 

-76 

50.88930 

450 

25 

33.38916 

- 76 

50.88012 

4^1 

25 

33.50436 

• 76 

50.87796 

452 

• 25 

33.60924 

• 76 

50.88396 

453 

25 

33.71250 

• 76 

50.87742 

434 

25 

33.80880 

^'76 

50.85204 

453 

25 

3'^. 37762 

-76 

50.82420 

457 

25 

33.94776 

• If 

50.71770 

458 

25 

93798 

- 76 

50. =3970 

459 

25 

33.92760 

'76 

50.55666 

500 

25 

13.90132 

-76 

50.45478 

-01 

25 

33.83208 

•76 

90.33406 

503 

25 

33.70278 

• 7 6 

50.11 098 

50U 

?5 

.’.3.69600 

-76 

50.03742 

500 

25 

33.70818 

- 76 

49.97136 

50f. 

25 

31.75396 

'-76 

49 . 9 39gf, 

5 17 

25 

’3.81 900 

.76 

49.91838 


2 5 

i.3. 92568 

- 7 6 

49.92390 

50- 

2^ 

34.36776 


45.94262 

5" 1 


.14.35056 

■ 76 

50.00904 

5 1 ? 


14 .48948 

- "^6 

50.04420 

s 1 3 

2 r 

■4.6 ?i<;? 

. 

50.08752 

51 4 

?✓ 

■’ 4 . 75939 

- 7^ 

50 . 11452 

1 5 


< 4 . s 6 7 4 y 


50. 1 0948 

S 1 

?':■ 

34.4462? 


50.10199 

SI 7 


I” . 301 2u 

n ' 

• I'l 

50 . .17558 

5 1 


■ 5 . t 14=2 

• 7 H 

49.96566 

5 


■' ^.,184 70 


i'9.8'7524 

-t 


■‘..0410 9 . 

■ ’’6 

49.73290 

> >> .^ 

- ^ 

34.94518 

1 76 

■J • 6 6 6 ? 

5/r3 

' “ » 

1 '.' .4 = 066 

‘ '’6 

sm 16 


2 ''' 


"76 

4<?e U»^h 6 ? 

•) >- 

PS 

. '•VU4 4 

■?6 


t •>, 


3 4 ,401 St, 

. 7 ^ 1 

1 r>5l u 

S -7 

'j. 

‘4.97-’2f3 

.7^ 


5 ^ 

'7 > 

.77424 

•- " 6 

iK. ^642i1 

^ t 


>.i . 8 ’ 860 

.t 7 k 

iA . 3f-l 


Table 1 coat’d 


SJPI 

531 

532 

533 

535 

536 

537 
535 

539 

540 

541 

542 

543 

544 
54C 

546 

547 

54ft 

549 

550 
552 
533 

554 

555 
536 
557 
555 
559 
600 
591 
602 

603 

604 

605 
6'J6 
607 

60ft 

604 

f 1 0 
61 1 
612 
s'3 
4 

1 

S’ 6 
617 
61.S 

Si 4 

j, 1 
•,'?? 
^ 


25 34.ft7o50 
25 34.96104 
75 35.06332 
25 35.19762 
25 35.47332 
25 35.58744 
25 35.70096 
2b 35.52420 
25 35.968ftb 
25 36.09444 

25 36.21222 
25 36.33456 
25 36.45222 
25 36.55993 
25 36.64902 
25 36.68994 
25 36.72930 
25 36.69216 
25 36.66672 
25 36.6427ft 
25 36.48444 
25 •«6. 38634 
25 36.28764 
25 36.16ftS4 
25 36.10110 

25 36.00372 

26 35.38342 

25 35.79102 
25 ■‘5.70670 
25 36.62920 
25 35.56230 
25 35.5129? 
25 -^5. 4733ft 
25 35.46ft22 
25 35.48820 
25 36.51274 
25 35.55270 
25 ^5.56636 
25 36.57214 
25 36.55995 
?5 36,52612 
25 36.47296 
23 35.39533 
25 36,30558 
25 36.22H64 
25 31.1340? 
25 35.03ftJ2 
■55 34.94136 
2'3 ‘4,34154 

25 •',4,75587 

25 34.69812 
25 34. 61 '24 
?s 34 ,54524 


. •‘‘J 








■^4 .4 I 

♦ 5 ' 


« 34U 24 




-76 48.86772 
-76 48.85716 
-76 48.85356 
-76 48.871°2 
-76 48.93516 
-76 48*95670 
.76 48.97326 


.76 ‘ 

48. 

99738 

.76 

49. 

03344 

-76 

49. 

05948 

-76 

49. 

07736 

• 76 

49« 

09650 

-76 

49. 

1 1 1144 

-76 

49. 

1 2482 

-76 

49. 

10946 

-76 

49. 

,06386 

.76 

49« 

,00578 

.76 

48- 

,Q0228 

.76 

48i 

.80322 

-76 

48 

.69858 

-76 

48 

.43566 

.76 

48 

.30318 

.76 

48 

.16572 

-76 

48 

.02634 

-76 

47 

.q2404 

.76 

47 

.79666 

-76 

47 

.66286 

-76 

47 

.53926 

.76 

47 

.41878 

-76 

, 47 

'.3164 = 


.76 47.201 10 

■,7ft 47.0833? 
-76 46.97304 
.76 46.871 10 

. If , 46.77486 
-.76 46.679?8 
.76 46.533’6 
.76 46.4933? 
.76 46*40448 
.76 46.32756 
-76 46.74800 
',76 46.19778 
-76 46.1 6304 
-76 46. 14198 
-76 46.1410? 

-7S 46.14996 

,7^ 46.16712 

'.76 46.17918 
•76 46.18446 
-76 46*19297 

■-7^, 46«i97'16 

-76 46*198‘'7 
-■’6 46.207''' 
-7^, 46.21320 
.76 46.25126 
-76 46.23612 
.76 46.24416 
^76 46 • ?6720 
76 46.26960 
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Table 1 coat'd 


f 


! 


■',' 1 . 2^432 
6-34 25 7 , 4.28958 

635 25 34.27338 

638 25 74.23760 

637 25 74.29384 

638 25 34.28503 

639 25 34.30674 

641 25 34.35126 

642 25 34.36032 

644 25 34.38840 

646 25 34.41864 

648 25 34.44672 


•76 46.26073 
-76 46.24134 
•76 46. .21878 
-76 46.22130 
•■76 46.21693 
-76 46.19880 
-76 46.19262 
-76 46.21086 
-76 46.23156 
-76 46.26228 
-76 46.29240 
-76 46.32318 
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Table 2 


TIME 

LATITLiDt 

1C135 

25 

34,16556 

1036 

25 

34.27326 

1037 

25 

34.38096 

1 w38 

25 

34.48866 

1 039 

25 

34.60068 

1 040 

25 

34.712'70 

1041 

25 

34.82634 

1051 

25 

35.997"^2 

1052 

25 

36.09744 

1053 

25 

36.18480 

1 054 

25 

36.26664 

1055 

25 

36.31536 

1036 

25 

36.33546 

1057 

25 

36.32844 

1058 

25 

36.29178 

1059 

25 

36.22998 

1 1 00 

25 

36.15084 

1101 

25 

36.05454 

1102 

25 

35.94300 

1103 

25 

35.82336 

1104 

25 

35.6=874 

1105 

25 

35.57484 

1106 

25 

35.45064 

1107 

25 

35.321 16 

1108 

25 

35.19j42 

1109 

25 

35.071 80 

1110 

25 

34.95438 

1111 

25 

34.84062 

1112 

73 

34.73706 

1113 

25 

34.64904 

1114 

25 

34.57494 

1115 

25 

34.51926 

1116 

25 

34.47504 

1117 

25 

34.44S88 

1118 

25 

34 . 431 60 

1119 

25 

34.421 28 

1 1 20 

25 

34.411 ^0 

1 1 21 

?3 

34,40904 

1 1 2? 

2^ 

34.40308 

1 1 ?3 

?3 

34.4152? 

1 1 24 

25 

34.42818 

1 1 2S 

25 

34,44966 

112^ 

25 

34 , 4791 2 

1127 

'>=, 

3u.,5'0852 

1 1 10 

?5 

;<4. 59.184 


tONRlTUOE 

^76 U6.'506>32 
-76 U.6. 88232 
-76 46.85832 
-76 46.83432 
-76 46.80414 
-76 46.77594 
-76 46.74546 
-76 46.39968 
-76 46.33324 
-76 46.36998 
-76 46.35372 
-76 46.36230 
•lb 46.37202 
-76 46.38510 
-76 46.40140 
-76 46.41960 
-76 46.44114 
-76 46.46202 
-76 46.48068 
'•lb 46.50348 
-76 46.53066 
-76 46.5624'’ 
^76 46.59726 
-76 46.63176 
-76 46.67124 
-76 46.71766 
-76 46.76520 
-76 46.80924 
-76 46.45112 
-76 46.89264 
'•lb 46.931 22 
-76 46.96860 

-76 46.99794 
.76 47.01924 
-76 47.03922 
-76 47.05254 
.-’6 47.05914 
.76 47.0626C 
^76 47.0596? 
.76 47.057=4 

:?6 :i7. 05980 
-76 47.06454 
-76 47.07225 
-76 47.08002 
.76 47.10324 
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Table 3 


TIME 

LAT i TUDE 

longitude 

319 

25 

43.50892 

-78 

1 1 .68164 

320 

25 

48.34800 

-78 

11.71710 

321 

25 

48.13702 

-78 

1 1.75256 

322 

25 

U8. 02610 

-78 

1 1 .73802 

323 

25 

47.8641 0 

-78 

1 1 .82222 

324 

25 

47.69850 

-78 

1 1 .85372 

32s» 

25 

47.52'’68 

-78 

11.881 50 

327 

25 

47.17698 

-78 

1 1 .92578 

328 

25 

47.00256 

-78 

1 1 .94354 

329 

25 

4(. .82724 

-78 

1 1 .95884 

330 

25 

46.65108 

-78 

1 1 .96952 

331 

25 

46.47468 

-78 

1 1 .97342 

332 

25 

46.31406 

-78 

1 1 .981 1 0 

333 

25 

46. I ^304 

-78 

1 1 .99040 

334 

25 

I! A. 01 443 

-78 

1 1 .99970 

335 

25 

4 -86556 

-78 

12.00456 

336 

25 

4 -.71 °64 

-78 

12.00696 

337 

25 

45.56994 

-78 

12.00810 

338 

25 

45 .42054 

-78 

12.00648 

339 

25 

45.24192 

-78 

1 1 .98890 

340 

25 

45.05526 

-78 

1 1 .96496 

341 

25 

44.86794 

'.78 

1 1 .93910 

342 

25 

44.681 88 

-78 

11.91402 

343 

25 

44.49378 

-78 

1 1 .89020 

344 

25 

44.30934 

-78 

1 1 .86764 

345 

25 

44.12358 

-78 

1 1 .8491 0 

346 

25 

43.9'461 2 

-78 

1 1 .84364 

347 

25 

43.79226 

-^78 

1 1 .84394 

348 

25 

4 3 .62600 

- n 

1 1 .84784 

349 

25 

43.45573 

W8 

1 1 .85285 

350 

25 

43.28904 

-78 

1 1 .86038 

352 

25 

42.96360 

- <'8 

1 1 .88630 

353 

25 

42.30424 

-78 

1 1 .90202 

354 

25 

42 . 63936 

-78 

1 1 .91624 

335 

25 

4?.46<426 

-78 

1 1 .93046 

356 

?5 

^^.303^4 

-7 1 

1 1 .94834 

358 

25 

4 > .97162 

-7f 

1 1 .98848 

35Q 

2S 

4 1 .80350 

-7d 

12.01 206 

400 

25 

41 .62674 

-78 

12.03216 

401 

25 

41 .45256 

,78 

12.05502 

402 

25 

4 1.2831.’ 

-78 

12.0832? 

403 

2^ 

ui . 1 Tiau 

-78 

12.11364 

404 

25 

4(1.9Aa56 

.78 

12.146/6 

4 .'5 


4M.3201 6 


11'. 1871 • 

■106 

25 

40.676‘»0 

•‘7-i 

1 2.228'. 

407 

25 

49 , S3So4 

^ 78 

12.27360 


2d 

40 . 3'^1’60 

-^78 

1 2.31758 

■1 VI 

25 

40.24314 • 

.78 

12.36006 

•1 1 11 

25 

4 4,00038 

•78 

1 2.40L2? 

4 1 1 

25 

5-4.0,' .1 :« 

-78 

12.44436 


.? ? 

> 1 , . 1 1 , 

-78 

12. 7766 

413 

25 

'.0 ,^7Ain 

.78 

12.00093 

I 4 

niA 

■> 9, .192 < 6 

^78 

1 2.54^06 

1 1 > 

25 

50 . 221)44 

• 7h 

1 2.88218 

4 1 

25 

<9. .46 

.78 

1 ?37H 

'll / 

25 

38,x84?6 


1 2 . 6 o 3 o 6 
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Table 3 cont ' d 


2*? 4 4. 

" > 4 4 
-44 

2v 4^ 


•7?37s 

.‘H 35s 
.25^7% 
.0921U 
.92090 
.75116 
.58090 
.U£,91 0 
.25358 
.11996 
.00812 
,.93678 
,.89^68 
,.901 86 
,.962‘*o 
'.04826 

49f -5 
' >2?« 
t. .1564 
<.01534 
<.31438 
< . 60580 
1.80R1? 
}. 18924 
1.47754 

1.77052 
'1.38450 
I. 40088 
'1.721 76 
I .34702 
I , 3?4i4 
1 .7 2740 
7 . 364 3 2 
v.40-,90 
2.11 63? 
3.2742? 
3 . 63256 

3.''<86'4» 
4.34566 
4 . 6 i 9 5 6 
4 191 ail 

5.3239-< 
5 . 655Tn 


-7^ 

1 ?.7 36a? 

-78 

12.75450 

-78 

1 2. 5 J256 

-78 

1 2.«4S70 

-78 

12.8902? 

-78 

1 2.92508 

-7ft 

12.96528 

• 78 

13.03380 

-78 

13.04298 

•78 

1 3.78816 

-78 

13.13796 

-7“^ 

13.19190 

-7b 

13.24950 

-78 

13.31364 

-78 

13.38300 

-78 

1 3.4530« 

-78 

1 3.52052 

-7ft 

1 3.58772 

-78 

1 *3.72944 

-78 

1 3.87524 

-78 

14.01720 

-78 

14.1 5340 

•78 

14.28672 

-78, 

14.41378 

-78 

1 4.54718 

-78 

14.67030 

-78 

1 4 . 7'9060 

-78 

14.91144 

^ 7ft 

1 S.0 35a<-' 

- 7ft 

15.16356 

-?ft 

! 5 . ~ 94 96 

. 7 ^ 

1 5.430 

73 

15.5691 

-78 

15.71 232 

' 78 

1 5.^5T6‘ 

-78 

16.37922 

^78 

1 ,4.93306 

7ft 

1 6. 3>'>a63 

-7.^? 

I 6.5440a 

-78 

16.70682 

-78 

1 6 . 8,4434 

^78. 

16.9921 ? 

*7.ft 

17.1 3984 

-78 

1 7 . 9.87p2 
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T»ble 4 


•^^5 

427 

428 

429 

430 

432 

433 

435 

436 

437 

438 

439 

440 

441 

442 

443 

444 

445 

446 

447 
443 

449 

450 

451 

452 

453 

454 

455 

456 

457 

458 

459 

500 

501 

502 

503 

504 

505 

506 

507 

508 

5?? 

Sts 

Si V 
s I ^ 
si ^ 


-AriTuoe 


|5 50.94618 
50.53696 
“ 50.40732 
25 50.22774 
:f? 50.06460 
49.73688 
49.57668 
27866 

If 49.12452 

If 48.97644 
If 48.32764 
If 48.65962 

If 48.49170 

If 48.32466 

If 48.14172 
If 47.97823 

If 47.81105 

If 47.63664 
If 47.46288 

If 47.28180 
If 47.09910 
-f 46.92570 
If 46.73280 
If 46.53420 
46.34940 
^f 46.1307 k 
I f 46.02340 
;f 45.8 o772 

If “5.71010 

I? 45.56172 

If 45.4297s 

••f 43.29196 

If 45.135^5, 

-f 44.97003 
25 44.80053 
25 44.6274.1 
If 44.4S040 

44.2695^ 

25 4«. 0723k 

If 4 '^.8781 8 

2p 4 .'i.6s« 64 
2r 43.4992? 

|5 43.31958 

;? 43.1 1 s?6 
If 42.76410 

4?.6?41 j 
;; 42.50s3;_ 

i. 4 -’.a (86s 

4?. 36216 

1^ 42.3631, 

_^f 4.'.40'17 k 
f" 42. 4688 k 

'42. 5694k 
4;>,^S244 
'■ 42.S209W 


long I tuoe 


.'78 
-78 
.-78 
-78 
-78 
-7p 
.-78 
'78 
.'78 
.'78 
-78 
.'78 
'78 
•78 
-78 
'78 
*7S 
.-78 
-78 
•78 
r78 
7 78 
'78 
.'78 
'78 
•78 
.'78 
.'78 
.'78 
-78 
-78 
'78 
'78 
'78 
-7« 

'78 
-78 

'78 
•78 
•78 
'78 
-7s 
•78 
•78 
.'78 

.-78 
-78 
-7,8 
-78 
-78 
•78 
•78 
•7s 
■78 
78 
78 


2.492<( 
2-6050 
2-6614 
2-7177, 

2.7815J 
2 - 8954£ 

2- 95080 
^•06744 

3*12048 

3- 17352 
3-22542 
3*26568 
3-30594 
3-34524 
^•37350 
3*41346 
3-45066 
3*47850 
3-50544 
3*52312 
3-55044 
3-57804 
3-59166 
3-60144 
3-62202 
3-65292 
3-69264 
3-72546 

3.76056 

3-80166 
3-85374 
3--70276 
3-93894 

3.97074 
4.00794 
4.03302 
4*06206 
4*08762 
4*10358 
4*12464 
4*14618 
4*16772 
4. 1 944? 
4-20960 
4.22790 
4.96348 
4*30350 
4*35066 
4.40244 
4-46094 
4-53552 
4-61454 
4-69572 
4*79326 

**.87,l5f, 

4.96572 
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I 


I 


I 


I 


r 


f 





f I ! I 


525 

526 

527 
‘=^26 

529 

530 

531 

532 

533 

534 

535 

536 

537 
535 

539 

540 

541 

542 

543 

544 

545 

546 

547 
543 

549 

550 

551 

552 

553 

554 

555 

556 

557 
553 
559 
60^a 
601 
602 

603 

604 

605 

606 
607 
:>03 
60Q 
^ 1 0 
6’ 1 

613 

I 'I 

n I y 
^ I o 
SI / 

->1 ^ 

M 

-.?l 


Table 4 cont'd 


25 

4? 

'.99=i04 

25 

43 

.17696 

25 

43 

.35035 

25 

43 

.51422 

25 

43 

.67112 

25 

43 

.82874 

25 

43 

.98126 

25 

44 

.11248 

25 

44 

.22786 

25 

44 

.35572 

25 

44 

.49258 

25 

44 

.63430 

25 

44. 

.78100 

25 

44. 

.92584 

25 

45. 

.07093 

25 

45. 

.22014 

25 

45 . 

.36048 

25 

45, 

.4912? 

25 

45. 

.6183b 

25 

45, 

,73470 

25 

45. 

,84732 

25 

45. 

96360 

25 

46. 

07Q40 

25 

46. 

20246 

25 

U6. 

33512 

25 

46. 

4»>952 

25 

46. 

61616 

25 

46. 

77066 

25 

46. 

‘53? 98 

25 

47. 

09736 

25 

47, 

26374 

25 

47 • 

41494 


25 47.56320 
25 47.70?40 
?5 u7. 83800 
25 47.^5123 
25 48.04098 
25 43.12882 
25 1+8.10830 
25 43.2395:> 
25 48.25080 
25 48.2131? 
?5 48.1440P! 
25 48.05010 
?5 47.93136 
2S 47.77914 
23 4 ' 7 , 61^78 
?5 47.44548 
?5 7^ 3r 

23 a 7 . 1 6 ? . 

i?3 46.9 1 558 
?‘> -I A. 7 n 36 
25 lih.3?31 
25 i6.j5?^3,1 

23 **^.1?>^48 

25 -J 5 . 3 0 M 
->5 'I 5 , 7^s 1 7 
'*5 U'.*)94,94 


-78 

5.07510 

-78 

5.1 8790 

-78 

5.29674 

-78 

5.40246 

-78 

5.50536 

-78 

5.60808 

-78 

5.70678 

-78 

5.79474 

-78 

5.87370 

-78 

5.95794 

-78 

6.04542 

-78 

6. , 3398 

-78 

6.22440 

-78 

6.31446 

-78 

6.40596 

-78 

6.49378 

-78 

6.58764 

-78 

6.67128 

-78 

6.75348 

-78 

6.83076 

-78 

6.90444 

-78 

6.97800 

-78 

7.05318 

-78 

7.1 3484 

-78 

7.22250 

-78 

7.31046 

-78 

7.40526 

-78 

7.50570 

-78 

7.61004 


7.71594 

-78 

7.81668 

-78 

7.90896 

-78 

8.00124 

-78 

8.08812 

-78 

ft. 1 75«4 

-78 

8.25600 

-78 

ft. 33724 


-7S 

ft. 42070 

-7H 

8. SI 148 

-78 

8.89902 

-78 

8.68812 

-7ft 

8.7646? 

-7ft 

ft . ft 3 7 1 6 

-78 

8.9041 8 

‘7ft 

8.96424 

-7ft 

9.01266 

-78 

9.0610ft 

-7P 

9.10416 

-7ft 

9. 1 r>024 

-78 

9.1 921 ? 

-7ft 
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Tcchnicnl Memorandum RS0 131 
Whiting Investigations on the 
Little Hahamas Bank 

Johr^ M. Hill*, Fabian C. Polcyn**, and Charles \’erri i 1 1 i on*** 

INTROmiCTinN 

•\n investigation ^'n the Bahamas to evaluate the 
capability of the high gain mode of Landsat 1 and 2 to deter- 
mine bathemetric data included an investigation of an anomaly 
called Whitings. The Whiting Investigation was condiicteil hv 
NASA/Coddard Space Flight Center, the environmental Research 
Institute of Michigan, the Remote Sensing Center c, .exa- 
Univeisity, the Applied Physics Laboratory of lohns Hopkins 
University, the U.S. Coast Cuard, and the ('ousteaia Society. 

The whitings were oi. served in satellite images 
even before *hc ship sailed. The whitings were investigated 
because they could conceivably cause an alarm due to their 
similarity to shoal areas. 

R/\CK"'.ROUNn ON POSSTBLF: ORICINS OF WIIITINCS 

Shn 1 1 OK- water sedimentation is a complicated pro 
cess. It involves chemical, plnsical and biological 

* Remote Sensing Center, Texas Afl'l University, College 
Station, Texas, 77 84S, l).S.\. 

** Fnvi ronmental Research Institute of Michigan, P.H. 

Box bl8, Ann Arbor, Michigan, 48107. U.S. A. 

*** MASA/C.odda rd Space Fliglit (Center, Creiuibelt, M.arvland, 

2 07 7 1, II. S. A. 
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relationships . The solubility of CaCO^ changes with varying 
levels in photosynthesis, respiration, precipitation of thitiO,, 
evaporation, rainfall and fresh water run-off from the shore- 
line (Bathurst, 11975J. Organisms not only secrete caibon 
skeletons; they are also a prime factor in bieaking and 
stirring them up. Tidal and wind-driven currents as well 
as wave action are the main factors in the abrasion and 
dissemination of these particles. 

The Little Bahama Bank, adjacent to the island of 
Grand Bahama, is similar to the Grand Bahama Bank in that 
it is a large flat bank surrounded and sheltered by rocky 
shoals. This creates an area of relatively undisturbed 
water. Water movement is important in the supplying of 
nutrients to biological communities and has an inf'uence 
on the distribution of carbonate which can be in solution 
(Smith, 1940; Costin, 1965; B. Katz, 1965, Traganza, 1967). 
Bathurst (1975) explains that tidal movements can be easily 
modified b> 'winds. From March or April to the end of August 
the prevailing winds are generally from the east or southeast 
and in September they begin to move toward the north. 

Surface water temperatures from 22°G to 31°(' 
(February-August) are attributed to warm waters brought to 
the Bahamas by the Gulf Stream (Smith, 1940; Cloud, 1962; 
Broecker and Takahashi, 1968) and to insulation on the 
shallow banks. <^i'-ud (1962) found that the lateral 
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variation on the Great Bahama Bank was only about ().5"C. 

The water in these shallow areas is relatively well mixed 
by the wind. 

Whitings consist of suspended aragonite muds 
that appear as clouds of milk-white water. The origii' ;f 
aragonite muds which make up whitings is currently under 
question and examination. Some in'^estigators believe 
these muds to be accumulations of aragonite needles from 
decomposed codiacean and dasycladacean algae ('.owenstam 
and Epstein, 1975; Stockman, Ginsburg and Shinn, 1967; 

Neumann and Land, 1969). 

Cloud (1962) conducted an investigation on the 
environment of calcium carbonate ceposition west of Andros 
Island, Bahamas which included a study of whitings. Me 
states that the propeller and dragging anchor of his ship, 
the P hysal ia , made whitings as they stirred up the bottom. 
Underwater springs as well as outbursts of gas can cause 
similar conaitions. Cloud also states that unusual local 
meteorological conditions may also produce identical effects. 
Ginsburg (1956) observed similar cloe.ds of sediment which 
were stirred up by schools of bottom feeding fish, such as 
mullet. These plumes arc termed "fish muds" and were 
observed in Florida Bay. Wells ai’d Illing (1964) investi- 
gated whitings in the Persian Gulf which they believe to 
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be related or similar to those found in the Bahamas. Tliey 
observed that the milky patches, or whitings, appeared and 
grew within a few minutes and persisted in the area for 
several hours. The visibility in these whitings was from 
less than 0.5 m to a few centimeters, the suspension con- 
sisted of aragonite needles and some pelagic organisms, and 
the density was approximately Ig of solid/lOOL. The 
following, although these observations sound as if the tail 
of a whiting instead of the source had been observed, are 
some interesting points made by Wells and Tiling: 

(1) The milkiness appears at the same moment over an 
area of many square kilometres. The region is 
normally devoid of fish which, even where they 
disturb the bottom, only muddy the water in 
patches a few hundred meters across. 

(2) The whitings form at the surface and overlie 
clear water. 

(3) They occur over any depth of water in the Gulf 
but arc very fr' .enc, though smaller, over dc|)ths 
of less than If* m. 

(4) The aragonite crystals flocculate, sink and disappear 
from sight in Jess than a day, whereas stirred bottom 
mud takes several days to sct’le. 
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(5) The underlying bottoir mud contains up to 
calcite, but in the whiti g there is "very 
little calcite". 

(6) Analysis of water immediately before and after 
the appearance of a whiting suggests that there 
is a very sliglit drop in Ca^^ and Ca^^/Mg^^. 

This could mean that calcium had been removed 
from the solution. 

(7) Water stratification in terms of temperature, 
salinity and pH remains stable and thi.s points 
to the absence of turbulence and bottom dis- 
turbance . 

Bahamian aragonite mud can be formed by either 
inorganic cr organic (physiological) processes. Stockman 
et al. (1967) have found that the algae are adequate sources 
of recent aragonite mud. A fact that is confusing is that 
there exist areas cf sheltered waters such as in tlic Bimini 
Lagoon that have prolific algal blooms with no aragonite 
muds (Bathurst, 1975). 

Cloud (1962) estimated that 5 wt. « of the mud 
is dctrital, 17-20 wt. I is skeletal i\':>ains, of wk.ich 4-5" 
is algal aragonite, therefore leaving 75 wt. " to be accounted 
for by some other process. This has led to the possibility 
of inorganic precipitation. Black (1 933) concludf'd that in 
shoal areas west of Andros Island the high rate of evaporation 
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generates a high concentration o^ salts, when in co.nbinatlon 
with the loss of CO^, causes precipitation of tiny aragonite 
crystal«:. Smith (1940) and Cloud (1902) showed that rates 
of withdrawal of CaCO^ were related to an increase of 
resiuciice time of the water on the banks, as well a'^ to the 
salinity and the accompanying fall in the produced aCa" • a(;('^ 
as aragonite is precipitated. Broeckcr and Takahashi (1906; 
state that the rate of CaCO^ precipitation is proportional 
to the degree of supersatui ion. T' Is data would tend to 
justify the inorganic origin of the aragonite muds, for this 
relationship does not neces‘'-ari ly have to exist if the proces 
is biological. They, however, could not account for a iciii^c 
percentage of the CO 2 budget which is lost from the Bank. 

It is assumed that he CO 2 ru.st be removed fro.Ti the Bank and 
is lost to the deep ocean. 

.\nother hypothesis (Weyl, 1961) is that in sea 
water saturated for aragonite, a sudden bloom of diatoms 
can remove so much CO from tne water that the only way 
to reach an equilibrium is by a widtsprf'ad i)rcci p itat ion 
of CaCO., thus causing whitings. DeGroot (1965), however, 
conducted laboratory experiments which tend to disprove 
Wells' and Illing's hypothesis. His findings indicated 
that sudden appearances of whitii.gs in the Persian Gulf 
could not have occurred because the mi mum rate of pre- 
cipitation of aragoni'c frcu sea water is mucli too slow. 
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l;vcfi fchcn hurricJ up in the laboratory, it still took at 

least tu 0 weeks to fore a precipitate. 

Aragonite iieeJle precipitates occurring at 

intervals of a few years also occur in the surface waters 

of the DcaJ Sea. See\ and Haery (19t»?) investigated one 

in which the suspended HCO. concent rat irn< first occurred 

at »he surface and progressed downward to a aaxiiaun depth 

of 40 n. T!»e hot torn sediments were found to consist of the 

18 1 > 

heavier isotopes Q and C '. The precipitation of HCO. at 

3 

the near surface, duo to evaporation, would cause the selective 
loss of lighter oxygen and carbon isotopes. Their data also 
favors the hypothesis of aragonite precipitation, but one 
should keep in mind that the water in the Dead Sea is not 
normal sea water. The extremely high salinities, which 
far exceed those found in either the Bahamas or Persian 
Gulf, could very well cause j'rec i p i tat ion of MGOj in the 
surface layers. Revel le and l-airbridge (1957) suggested 
that temporary morning clouding in still waters in the 
lagoon of lloutman's Ahrolhs /southwestern Australia) might 
be due to photosynthetic reduction of Cu^ pressure. 

(,’loud (19().’) studied three whitings just west of 
Andros Island. The first whiting drifted away before it 
could be entered, but an underwater inspection of the 
drifting margin proved the water to be unusually milky due 
to suspended matter. No local turbulences or bottom dis- 
turbance.s were noted at the site. 
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Ur>on entering the second whiting. Cloud states 
that ’’its center gave the sensation of weightless fixity 
in the nidJle of a sunlit cloud hard. it was iapossihlc 
without resting notionless, to detect buoyancy and drift 
or to tell up from down or sideways. The brilliant lighting 
was so dispersed that a hand, invisible at arn's length, 
had to be extended to grope for bottom and avoid a collision 
on surfacing. Xo fish were seen (or felt) nor was any other 
evidence of bottom disturbance found". The suspended material, 
upon microscopic examination, was found to consist of primarily 
aragonite needles. 

The third whiting was about 4 km long and 0.8 km 
wide. During a six-hour study no large schools of fish were 
observed except for a few small Hal ichoere s , tiny Eques found 
in and about the sponges and much later five sharks were seen. 
No other whitings were seen to occur in the surrounding area. 
Chemical analyses of water samples from inside and outside 
the whitings were examined. The whiting water, after particles 
larger than 0.45 microns were filtered out, was significantly 
high in calcium, phosphates, CO 2 evasion, pH, the rate of 
photosynthesis, and low in alkalinity and partia-’ piessurc of 
CO^. The whiting kept its properties for at least 45 hours 
without being eliminated by diffusion or from mixing caused 
by external sources. 
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rioud noted that all of the whitincs were elonr.ate-.i 
and drifted with the wind and tidal current. school inc 

of fish was observed in the area, the whitincs did not drift 
away from a fixed point or move around erratically, and no 
unusual neteorolop ical or other disturbances were noted. 

He made no definite conclusions as to the specific cause 
or origin of these whitings except for the mentioning of 
the possibility of the previously stated ideas of inorganic 
or biological processes. 

Cloud stressed the need to examine individual 
whitings in any given area to determine what features 
they cither have in common or in which they differ. 

KMITING OBSRRV' ATIONS AS COLLF.CTl-n 
FROM THH CALYPSO 

Two whitings were investigated by the scientists 
and crew of the R/V Calypso on the fifth of September, 1975 
on the Little Bahama Bank. Whitings occurring in this area 
were observed in satellite images even before the cruise. 

Upon entering the area, a T-38 aircraft was utilized 
to obtain a synoptic view of the area and to locate a group 
of whitings. Once the area of whitings had been identified 
by the jet, the Calypso was directed to the source of the 
whiting by the Calypso helicopter. A Zodiac, a small rubber 
raft, was deployed from the ship and sent to the source. 
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Upon arriving at the source the Zodiac was anchored and 
the divers proceeded to enter the whiting. The Zodiac was 
used hecnusc of its <;nccd, maneuverability, and it lessened 
the chance that schools of fish would be scared away if 
present. The excellent naturalistic instincts and experience 
of Captain Cousteau’s divers made it possible for several 
new and interesting observations. 

The Zodiac crew collected measurements of depth, 
sccchi depth, sea surface temperature, surface salinity, 
and water samples for later analysis. The divers collected 
bottom temperature, bottom salinity, water samples, sediment 
samples (Figure 1) depth vertical visibility, horizontal 
visibility, current direction, a rough estimation of current 
speed and general bottom descriptions and comments (Table 1). 
Water and sediment samples and observations were collected 
from the source, the tail, and from the clear water surround- 
ing the whitings. 

The source of the first whiting. Station P.B, had 
several sharks swimming within it. One shark even took a 
swipe at a diver's flipper. Investigators have been known 
to not enter the whitings for just such a reason. Once in 
the plume the divers observed Cal 1 1 anassa mounds, with 
heights of 6 to 10 centimeters, covering the bottom (Figure 2). 
Bathurst (1975) states that conical mounds constructed by 
Call ianassa , burrowing shrimp, are common in stable sand 
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Fig. 1. Diver collecting sediment samples for later analysis. 
Note the very fine sediments. (Courtesy of 
Joe Thompson and the Cousteau Society) 










Fig. 2. Callianassa mounds surrounded by blades of Thalassia . 
(Courtesy of Joe Thompson and the Cousteau Society) 
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TABLF 1. Data Obtained from the Whiting Investigation on the 
Little Bahama Bank by the R/V Calypso (September 5, 


1975) 



F.B 

Source 

F.A 

Out of 
Whiting 

F.C 

Source 

ED 

Tail 

Depth (m) 

7.5 

7.0 

6.5 

7.0 

Seech i (m) 

7.0 

Seen on 
bottom 

3.75 

4.25 

Surf. Temp. (“C) 

29.93 

29.94 

30.98 

31.20 

Bot. Temp. (°C) 

29.90 

30.13 

30.24 

30.19 

Surf. Sal. (°/oo) 

38.294 

38.307 

38.457 

38.348 

Bot Sal. (°/oo) 

38.316 

38.319 

38.454 

38.342 

Vert. Vis. (m) 

2.0 

7.0 

1.0 

2.0 

Hor. Vis. (m) 

0.5 

5.0 

2.0 

2.0 

Cur. Dir. 

Minimal 

S.F. 

S.F.. 

S.E. 

Cur. Sp. (knts) 

Minimal 

0.25 

0.75 

2.00 

Mounds 

Present 

Present 

Present 

Present 

Sharks 

Present j 

1 

None 

Observed 

None 

Observed 

None 

Observed 

Fish 

Present 

None 

Observed 

None 1 

Observed 

None 

Observed 

Meteorol igical 
Disturbances 

None 

Observed 

None 

Observed 

None 

Observed 

None 

Observed 
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habitats. The usual mounds resemble miniature volcanic cones, 
have remarkably constant sizes, have a basal diameter of 
about 20 cm, a height of about 6 cm, and have approximately 
30“ sloping sides. At the peak of the mound there is a vent 
about 3 mm wide. The cone has a surface of loose sand. Hrain 
have occasionally been seen being expelled from the exhalant 
vent. Shinn (1968) originally discovered that the mounds 
were constructed by the crustacean Callianassa . The bases 
of the closely arranged mounds were surrounded by blades of 
Thalassia . 

Ginsburg (1975) explained that these shrimp 
homogenize the sand to a point where '‘the history of 
deposition and with it the story of successive sea-floor 
environments is hopelessly jumbled". Bathurst (1975) 
remarks that in some areas the mounds are the only source 
of loose grains. Tidal currents derived from evidence 
obtained around the Berry Islands can move the loose grains 
about 3 to 4 cm at the most and with time, that could be 
significant . 

A large school of what was termed Bone fish were 
seen bouncing off the bottom. These fish most probably had 
attracted the sharks into the area. The fish could have been 
either bottom feeding or merely cleaning parasites or other 
growths from their undersides. In this particular area of 
the whiting, the source was covered with a clear surface 
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layer approximately 2 m thick. The surface currents, usually 
beinfi more rapid than bottom currents, could have kept the 
sediments from rising to the surface at the source. The 
sediments in the tail of the whiting would eventually disperse 
into the surface layers. The water was turbid from this clear 
layer all the way to the bottom. 

A diver that had left the Calypso, which was 
anchored near the middle of the whiting, observed relatively 
large amounts of sediments being expelled from the Cal 1 i anassa 
mounds. This observation raises the question as to whether 
or not the Callianassa may be a possible cause of whitings. 

In the clear water surrounding the source of the 
same whiting. Station EA, the bottom was again covered with 
closely arranged Callianassa mounds surrounded by Thalassia 
blades. No sharks, fish school, or burrowing shrimp ivere 
seen disturbing the bottom sediments. New sources of whitings 
were observed by the crew of the helicopter. The new sources 
appeared as billowing, or mushrooming clouds of white milky 
sediments that were being forced to the surface. From search- 
ing the literature this is quite possibly the first actual 
observation of the evolution of a non man- induced whiting in 
the middle of a shallow bank. It is very difficult to 
imagine what mechanism might mushroom bottom sediments to 
the surface in about 7 m of water. That is, in fact, if 
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the sediments do actually come from the bottom. The large 
sharks which were observed in the source area of the first 
whiting could have created a turbulence in the water column. 
Rczak (1976) , in conjunction with the Shell Development Cor- 
poration, has studied sediments coJlected from whitings in the 
Persian Gulf with an electron microscope. His findings were 
that the sediments consisted of primarily aragonite needles 
and fine, fine fragments of mollusk shells which had to be 
stirred up from the bottom. Currents roughly measured in the 
whitings ranged from about 0.25 to 0.75 knots. At one station 
the current was negligible. These values are about the ex- 
pected for the area. The salinities both top and bottom 
ranged from a mean of 38.358“/oo to 38.352®/oo respectively. 
These summer values are rather high when compared to open 
ocean waters which are usually around 36“/oo, but values of 
46“/oo have been recorded just off Andros Island (Cloud, 1962). 
The temperature and salinity values point out that the water 
mass was fairly well mixed from top to bottom and that the 
possibility of underwater springs, that could have caused these 
particular whitings, was not probable. 

Station F.C was located at the source of the second 
whiting. This was the most dense whiting source because, the 
secchi depth was only 3.75 m. The bottom was again found 
to be covered with mounds and surrounded by Thalassia grass. 

Ill 








T 


I 

w 


J, 




.. . 



' 1 1 ! I ? I I 

« 


1 


I 

i 

f 

i 





i 



The sediments between the mounds appeared to be darker 
than the previous area. No sharks, schools of fish, or 
shrimp were observed stirring up the sediments. The 
absence of any of the previously observed sources of 
disturbance from the first whiting left speculation. 

No meteorological disturbances were observed in the local 
area over either of the two whitings investigated. 

Station ED was established to investigate the 
tail of the second whiting. The sediments ivere still quite 
concentrated even when about 1 km from the source since 
the secchi depth was 4.25 m. The mounds surrounded by 
grass were again observed. No sharks, fish, or shrimp were 
seen in the tail of the plume of the second whiting. 

The general shape of these whitings can be con- 
structed from the observations (Figure 3). A layer of 
clear water was over the source of the first whiting. The 
sediments were also the most dense from top to bottom in the 
source of both whitings. The vertical visibility observa- 
tions made on the bottom, in the tails of both whitings, 
especially in the first whiting since it was most probably 
the oldest and longest of the two, was much better. The 
better bottom visibility and the lack of a clear surface 
layer would lead to the belief that the heavier sediments 
have already settled out and that the finer particles are 
carried to the surface layers in the tail of the whitings. 
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CONCLUSIONS 




Several points of interest were made by the 
investigation of these two whitings on the Little Bahama 
Bank . 

(1) The whitings were all elongated in form. 
(Figure 4). 

(2) The source areas of the whitings did not 
move significantly from a fixed point. 

(3) The suspended particles appeared to be 
evenly distributed through the vertical 
water column, except for a 2 m layer of 
clear water on the surface at the source. 

(4) They displayed definite source areas from 
which the sediment was dispersed in a nearly 
straight course by the wind and tidal current. 

(5) The whitings kept their discrete properties 
for long distances to the point where 
several whitings were observed to be parallel 
for relatively long distances while being 
separated by only 4 to 5 m or less of clear 
water (Figure 5). They would eventually mix. 
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Fig. 4. Elongated shape of the whitings and the R/V Calypso. 

(Courtesy of Joe Thompson and the Cousteau Society) 
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Parallel courses of several distinct and s'jparate 
whitings. (Courtesy of Joe Thompson and the 
Cousteau Society) 
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(6) Two possible new causes of whitings have 
been observed: 

(a) A large school of fish, accompanied 
by sharks, was observed stirring up 
the bottom sediments in the source of 
the first whiting. 

(b) Cal 1 ianassa , burrowing shrimp, mounds 
ivere observed dispersing sediments to a 
height of a meter or more in the middle 
of the second whiting. 

(7) The formation of a new whiting was observed 
aerially from a helicopter and appeared as 
billowing, mushrooming sediments being 
presumably forced from the bottom to the 
surface by an as yet unexplained source. 

Further sediment analyses being conducted by 
Ginsburg at the University of Miami will hopefully add to 
the knowledge obtained on the origin of these two whitings 
on the Little Bahama Bank, It is also hoped that these 
observations may generate future studies on the origin 
of whitings that will take into consideration the old as well 
as these new hypothesized sources of the aragonite particles. 
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GKNRRAI, SAMPLING INFORMATION 




Zodiac (Rubber Raft) Samples : 

(1) Seech i 

(2) Sea Surface Temperature (Bucket Tliermometer) 

(3) Surface Salinity Samples 

(4) Gallon Surface Water Samples 

Diver Collected Samples : 

(1) Bottom Temperature (Reversing Thermometer) 

(2) Bottom Salinity Sample 

(3) Gallon Bottom Water Samples (Niskin Bottles) 

(4) Bottom Sediment Samples 

(5) Depth 

(6) Vertical Visibility 

(7) Horizontal Visibility 

(8) Current Direction 

(9) Current Speed 

(10) General Bottom Description and Comments 
Helicopter and T-38 Jet : 

Spotted source and tail of whitings. Also spotted clear 
areas out of the plume. 




DATE: 9/5/75 

STATION: RA 



* The samples were collected from clear waters outside of 
the plume or whiting #1. 

Sample Code: 

(1) AT-1 Gallon Water Sample Collected From the Surface 

(2) ATS- Salinity Sample Collected From the Surface 

(3) AB-1 Gallon Water Sample Collected From the Bottom 

(4) ABS- Salinity Sample Collected From the Bottom 

(5) A,A,A-2 Sediment Sample From the Bottom 
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Secchi : Seen laying on the Bottom ■ 

Sea Surface Temperature: 29.94°C , 

Botcom Temperature: 30.21°C-(A = 33.2®C) (possible diver problem) 

Depth: 7 m 

Vertical Visibility: 7 m 

Horizontal Visibility: 5 m 

Current Direction: From SE , 

Current Speed: Material moved 98" in 19 seconds 

Comments: 

(1) Bottom covered by mounds 4-12" high 

(2) Eel crass surrounded mounds 

(3) Three bottom samples were collected from top 
of mounds 





(4) No sharks were observed 

(5) On bone fish were observed 

(6) No mollusks were observed 

(7) Helicopter pilot spotted new sources that 
appeared to be billowing (mushrooming) 
clouds of sediments coming to the surface. 




DATE: 9/5/75 

STATION: EB (Source of First Whiting) 



(1) BT-2 Gallon Samples Collected on the Surface 

(2) BTS- Salinity Sample Collected on the Surface 

(3) BB-2 Gallon Samples Collected From the Bottom 

(4) BBS- Salinity Sample Collected From the Bottom 

(5) BBB-3 Sediment Samples From the Bottom 
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Secchi: 7.0 m 

Sea Surface Temperature: 29.93°C 

Bottom Temperature: 29. 92°C - (A = 30.1®C) 

Depth : 7.5m 

Vertical Visibility: 2.0 meters above bottom 

Horizontal Visibility: 1.5 Feet 

Current Direction: No Current 

Current Speed: No Current 

Comments : 

(1) Mounds from 4 to 12 inches covered the bottom 

(2) Eel grass surrounded the mounds 

(3) Three bottom sediment samples came from the 
top of the mounds 

(4) Sharks were observed near the Zodiac 

(5) Bone fish in large numbers were observed 
hitting and bouncing off of the bottom 

(6) Mollousks were also observed by one diver to 
be squirting up sediments 

(7) The source was spotted by helicopter 

C8) There was a layer of relatively clear water 
(2 in thick) on the surface, but was turbid 
from 2 m to the bottom 
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D/VTE: 9/5/75 

STATION: EC fSource of Second New Whiting) 



Sample Code: 

(1) CT-1 Gallon water sample collected from the surface 

(2) CTS- Salinity sample collected from the surface 

(3) CB-1 Gallon water sample collected from the bottom 

(4) CBS- Salinity sample collected from the bottom 

(5) C,C,C-3 Sediment samples collected from the bottom 

(See Explanation) 

(6) CO-1 Sediment sample collect from the bottom outside 
of the plume 




Secchi: 3.75 m 


Sea Surface Temperature: 30.98°C 

Bottom Temperature: 30.3°C (A - 32.3°C) 

Depth: 6.5 m 

Vertical Visibility: 1 m before bottom was visible 

Horizontal Visibility: 2 m lateral visibility 

Current Direction: From S.E. 

Current Speed: 9 seconds to go 18 inches. 

Comments : 

(1) Bottom covered by mounds 4 - 12" in height. 

( 2 ) Eel grass surrounded mounds. 

(3) 2 bottom sediment samples collected from mound 

(4) 1 bottom sediment sample collected from between 
mounds 

(5) No sharks were observed 

(6) No bonefish were observed 

(7) No mollusks were observed 

(8) Sediments between mounds looked muddier and 
darker, mound sand looked clean and white. 
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RATE: 9/5/75 

STATION: HD (Tail of Second Plume) 



Sample Code: 

(1) DT-1 Gallon Water Sample Collected From the Surface 

(2) DTS- Surface Salinity Sample 

(3) DB-1 Gallon Water Sample Collected From the Bottom 

(4) DBS- Bottom Salinity Sample 

(5) Sediment Samples: 

D = From top of Mound 

D+ * From Between Mounds 






D++ = From Top and Bottom of flounds 


DO = From Outside of the Plume 
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Secchi: 4.25 m 
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Sea Surface Temperature: 31.2®C 
Bottom Temperature: 30.20 - (30.25-A) 

Depth: 7 m 

Vertical Visibility; 2 m from bottom 
Horizontal Visibility: 2 m 

Current Direction; From SE (45°) 

Comments : 

(1) Bottom covered by mounds 4-12" in height 

(2) Eel Grass surrounded mounds 

(3) Three bottom sediment samples collected 

(See Previous page for sample code) 

(4) No sharks were observed 

(5) No bonefish were observed 

(6) No mollusks were observed 
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DATE: 9/5/75 

STATION: RE 

One Bottom Sediment Sample Was Collected About 5 miles 
From the Sited Whiting Area. 

Sediment Sample Is Marked '' + ++". 
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